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ABSTEACT 


Linear Time Invariant digital filters are charact- 
erized by difference equations with constant coefficients. 
They generally require complex hardware structures when 
they are reaxised physically. This thesis deals with the 
development of a non-linear digital filter that is 
equivalent to a linear digital filter, in the sense that 
the expected values of their outputs are equal. It is 
shown that for a fixed input sampling rate, the power 
spectrum of the output of the nonlinear filter can he made 
■^o asymptotically approach the power spectrum of the 
output of the linear filter. The hardware requiriment 
for the nonlinear filter is simpler than that required for 
the linear filter, pcssihle application of this filter 
structure is in the realisation of filter hanks, and in 
the realisation of complex filters, such as a fading . 
dispersive channel simulater. 
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CHAITER t 
INTRODUCTION 

Digital signal processing techniques find increasing 
application over analog signal processing techniques owing 
to their greater accuracy, reproduciblity and programmability . 
/idvances in integrated circuit technology have made, for the 
same accuracy, dsp techniques more economical than analog 
si'cnal processing techniques. Tlie theory of dsp emerged 
from an understanding of sampling and its spectral effects 
and the use of the z-transform. The integrated circuit 
implementation of digital circuits caused dsp to further evolve a 
and it. continues to do so today. 

The sampling of a signal is quantization of the 
Signal in time. Similarly the concept of the quantization 
of the signal in amplitude developed. The late 40 's saw the 
introduction of the statistical theory of amplitude quanti- 
zation. One of the numerous uses of quantizers was in 
correlation devices. Digital signal processing used finite 
wordlength machines, which necessarily meant quantization 
in amplitude. The statistical theory of amplitude 
quantization was applied to study the finite wordlength 
effects in digital machines and, in dsp areas, specifically 


in digital filters. 



Linear digital filtering ia one area of digital 

signal processing, v/l.^-^re the digital filter is the 

iriplemontation of a constant coefficient linear difference 

equation. The hardware realisation of a digital filter 

therefoi-ft requires a number of basic hardware units e.g. 

multipliers , delay elements, and simmers, the number of 

elements increasing as the order of the filter, hence 

% 

increasing the cost. Kirlin [7j introduced the statisti<y>l 
coarse-sampled-data filter with a view to reduce the 
hardware requirements of the linear filter. It was his 
contention that the performance of the coarse-sampled- 
data filter matched that of the linear digital filter It 
meant to replace. 

Kirlin 's nonlinear filter appears to be very 
promising, for while it is said to be equal to the linear 
digital filter in its performance, its hardvxare requiremonta 
are simpler, impli/ing a reduction in hardware complexity 
as well as cost. This aspect is particularly important 
in the realisation of higher order systems. For e.g.» 
a digital hardware unit may simulate a fading dispersive 
channel for eg^ a troposcatter channel. The large bank 
of filters required to adequately simulate multipath delays, 
and the multiplication by complex coefficients to simulate 
finyleigh fading and Doppler Spread, require a large amount 



ot comilex hardware. Implementing a fading dispereive 
channel Simula tor using Kirlin*a scheme would mean a 
big saving in hardware. 

The object of the present work . was to evaluate 
the performance of Kirlin's coarse-sampled-data filter 
with respect to that of the linear digital filter. If the 
performance of the coarse-sampled-data filter was found to 
be not equal to that of the linear djLgital filter, could 
the structure of the coarse-sampled-data filter be chantged 
so that the performance of the modified filter matched 
closely, that of the linear digital filter? 

It was found that the performance of the coarse-samplod 

daf"' filter when the input was a vrhite noise sequence, was 

nuv-rhere near that of the linear digital filter. However, 

of 

altering the structure / the feedback of the coarse-sampled- 
data filter led to ^diat has been called 'N-bit feedback 
sampled -data filter'. The performance of the output of the 
N-bit feedback sampled -data filter, when the input was a 
white noise sequence, was found to asymptotically 
approach tliat of the linear digital filter. 

Organization of Chapters : 

In Chapter 2, theoiy and techniques of digital 
filtering, the statistical theory of quantization and the 


coarse-sampled-data filter using an RC averager, are 
reviewed. In Chapter 3, the performance of the first 
order and second order linear digital filters, coarse- 
sampled -data filters and N-bit feedback sampled-data 
filters are evaluated, for a white noise sequence input. 
In Cha.pter 4, the performance of first and second order 
Il-bit feedback sampled-data filters as obtained by- 
simulating the filters on the DBC-10 is discussed. 
Chapter 5 summarises the results obtained in Chapter 



Chapter 2 

REVIEW OF DIGITAL SIGlIAL TROCESdl HG AND T.IB STATISTICAL 

THEORY pF QUANTISATION 

In this chapter, digital filter design techniques, 
effects of quantization on analog signals, use of quantizers 
for autocorrelation determination, and the statistical coarse 
sampled wjata filter are discussed. 

2.1 RELATION BETV/EBN ANALOG AND 'DISCRETE SYSTEMS 

The Fourier tr'>nsform of an analog signal is given 

by 

•A 

. . 6-3“’^ d-t (2.1.1) 

Inversely is derived from by the relation 

i 

^ jx^(ato) e^^ d.*> (2.1.2) 

— 3C 

If be sampled, with sampling rate 1/T, it can be shown [l 

that for the resultant signal x(nT), the Fourier transform is 
given by 

Kjc*-') = f- ^ (2.1.3) 

k=-oo 

where is the sampling frequency in radians and oa = 2ti/T« 

O 0 

X(jo>) is obviously periodic with period 2tt/T. 
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Lf is band limited over the ran/-;<» {-n/T, n/T)^. then 

X^(j<«) r.n] X(j'o) are Identical in this frequency range. Other- 
wise a condition known as 'aliasing' occi.’ra. 

When is a stationary stochastic process, it can be 

shown [2], that the ensemble autocorrelation function R (1) 

A 

of the sequence x(nT) an! the ensemble autocorrelation function 
R (j) of the analog signal x. (t) are related by 

Xa a 




1 t=1T 


(?.i 


It can also be shown [2] that if is ergodic in the 

Autocorrelation sense then the time averaged autocorrelation 
function of XA(t), defined by 

£ 

=■ 2K^ J St 

A 

and the time averaged autocorrelation function E^(i) of the 
sampled signal 3?(nT) defined in like manner, are related ws 
follows : 


rIO) = RJ ( < ) 1 (2.1.5) 

The power spectrum and the ensemble AGP of form a 

Fourier transform pair thus : 


and 


^A 


r 


*A . 


R^ (x) dT 




J 


( 2 . 1 . 6 ) 



IThe power spec. " r j h of the sa^”- 




X 


Sjz) 




/ _ R, Mr) 55 


-k 


5&=e 


ji) T 


(?.u o 


.£[R^(m) T') 


P^qr. (2.1.4) cTfectively implies that 'i„(l) may be obtaint-i bv 
sampling A (t) at T second intervals and because 
R ( S ( 3<2>) 


oo 


. 3..(a^>:»T) =£[R_(k)] =i 2- B ( 2 . 1 . 8 ) 

A J. <90 A, 


by iuvol ing the relation (2.1.?). 


2.2 DIGITAL i-'ILTER DBSTGW TECHNIQUES 
2'. 2.1 Filter Properties ; 


Important properties of digital filtets are magnitude 
sq'sured response H(e''^) , phase delay T^OQ), group or 

envelope delay T (^^'') . 

O 

If the filter transfer fun'^tion is given by then 


Tp( ) 


[il(z) 


Tan 




UJ 


Ini(H(z)) 

Tmdrr 


z-e' 


] 




In 


rlilgi.. ! 
) 1 


] 


Ji 


(r P. 1 ) 

(2.^,2) 


(2.2.3) 
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Tp(jto) = 


( * 

z I ; 
I z=e‘ 


(2.2.4} 


, I Z *“0 


(2.2.5) 


= - fie[z (In H(z))] (2.2.6) 

Those poles of H(e^^') ^ that lie inside the unit circle in 
the z~plane are determined to be ‘‘'•c'poles of H(z). The zeros 
of H(z) are not uniquely determined, but choosing^ them to be 
tho'^e zeros of H(e"^^'' ) ^ that lie inside or on the unit circle, 
yields a miniroiam phase filter, dome si^nel processing applica- 
tions req.uire the filter to be linear phase i.e. that tiie group 
delay b6 zero or a constant tjver the frequency bands of interest. 

2.2.2 finite IXirafclon Impulse Response Filters : 

Tvro widely used filter types are FIR and ITR filters. 

Rcurip advantages of FIR i i.ltf-rs are that they are always stable 
and using an appropriate delay can always be reali/.ed. They 
can be designed to yield exactly linear phase, a (-ondition 
required for speech signal processing an'’ data transmission. 

Disad vantages of FIR filters are that a large value of impulse 
response duration is required to adequately approximate /diarply 
cut off filters, hence requixing a large amomit of computation 
and time. But with tha advent of the FFT, fast convolution has 
been made possible, so that FIR filters are now competitive 
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with sharp cut-off IIT filter designs. If an FIR filter with 
impulse response duration = N, be required to be linear phase 
with T (j/j..) = -a then 

Jr 

a = ^ ' ■ (2.2.7) 

and h(n) = h(N»1~n) (2.2.8) 

where jh(n), n = 0, N-l is the impulse response 

sequence. 

FIR Filter design techniques : 

(a) Windowing ; The Fourier Transform of the impulse response 
of a digital filter is periodic and it can be expanded in a 
Fourier series thus ; 


430 


H(ej“)= Z h(n) 


( 2 . 2 . 9 ) 


n=- oc 


where | h(n) J is the impulse response of the filter and 
given by 


2TC 

h(n) = ^ dfu , ~a.<n 

o 




( 2 . 2 . 10 ) 


Because ] h(n) t extends from -£>;>< n t the filter is 

I ; 

unrealizable. A finite weighting function w(m) is used to 
modify the Fourier coefficients | h(n)j' so that 
■^(n) = w(n).h(n). Then if H(e^^ ) h(n), H(e^^*'' ) should 
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aj'proxiiaate ) within some paraissible error. Clearly 

) is the circular convolution of ¥(e^'^^) eni 
.lfW(e^^ ) = 6(tb>) then H(e^ ) = K{e^ ), but if w(n) is a 
finite sequence, W(e^^) is not an impulse. Then ¥(6^*^*^ ) ■’iust 
be chosen such that its main lobe has a small width and hence 
hi;’;her energy concentration; the side lobes should decrease in 
energy rapidly as Co -- 7 > %, So:'‘^' windows use'? are Rf'-ctr^ngular 

windovr, Hamming window, Kaiser window; 


(b) Frequency Sampling : Given N equally spaced points 
in the TjFT of an impulse response sequence, the 
transfer function H(z) can be foxmd [l, pp 10^-107] to be 


liU) “ 


-N 


N-l 

■‘=^r 

k=0 


Hii£l 


( exp(-> ) 


( 2 . 2 . 10 ) 


Here the transfer function n(z) approximates a continuous 
transfer function. The approximation error is zero at the N 
specified frequencies and finite between them. 

(c) Optimal (Min Max Error) Filters ; Considering the 
linear phase FIR filter design problem as a Ch^byshev Approxima 
tion problem, a set of condition is derived for which it can 
be proved that the solution for which the peak approximation 
error over the interval of approximation is minimized , is 
Optimal. Dv^fine 3^(6*^ " ) = the transfer fn, being approxi- 

mated, '.5j(e^^'^^) = a weighting function on the epproxln-a tion 
error, so tiiat the, relative size of the error in diff<=*rent 
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*1 ^ 

frequency bands may be chosen, HCe** ) = the frequency 
response of the designed filter. Then 


E(e 


T 


) = 




tH„(e 




( 2 . 2 . 1 ; 


i T 

where ) is the approximation error* Then the Chebyshev 

3 pproximation problem may be stated as 




= min [max ] 

min coeffs ioeA 


( 2 . 2 . 1 ,^ , 


where A is the disjoint union of frequency bands of 1ntei 


Of these three design technique, the frequency sampling and 
Minmax error techniques are superior to windowing, 

2.2,5 HR Filter Design ; All filter design techniques shoul' 
leod to stable and realizabj.e filters, therefore HR filters 
'should he such that 


h(n) = 0, 


h(n) 


n < 0 


n) y o 


(2.2.t4 / 


The general form of the HR filter transfer function is 


M _ N 

h(z) = (21a^z~^)/(l + X- h^z”^) (2. 2. IS) 

i=0 i=1 ^ 

where at least one b^ is nonzero and the roots of the denomina^ 
tor are not cancelled by roots of the numerator, which case 
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leaif, to an FIR filter. For der’ive^! from an analog 

f ilter, 114 ; N. 

;’e3i/rn Techniques : 

Inipulae Inva.rian.t Transformation ; Using this transforraatton. 
tUe impulse response of the resulting digital filter i« a 
sampled version of tha.t of an analog filter. let H(6) have 
simple poles and M ^ N, then 

Sr- 0 * , 

' ''(s) = 21 ^ 1 0., = H(s) (ai-d.)/ (2-2.16) 

i=1 ®'^^i s=-d^ 

The analog filter imt'ulse response is given by 

IJ ^ 

h(t) = 7" *^4 e ^ . u (t) (2.2.17) 

■ • i=1 

■ I 

/. h(nT) = 7 e ^ u_^^(nT) (2.2.18) 

i=T 


whore T is the sampling period. 

N . i\ 

Eiz) = 2~ h(nT) 

i=1 i=1 l-e Si 


(2.2.19) 


~d.T 

Obviously, this uses the mapping, s+d^—^* 1~o z"" J for a 

i iO T 

simple pole at ~di. Since H(z=e‘' ) is the Fourier Transform 

of the sampled signal h(t). 

O 0 

H(e^‘^^) = i fL H( j (S'-^F’wT) , where w = 2 ti,/T. 

^ k--=~oo ® ® 



to an i'lR filter, i'or H(-?) derive'l from an analog 

f ! Iter, I'j F. 


Techniques : 


Impulse Invariant Transformation : Using this transforcss tion, 
tre impulse response of the resulting digital filter is a 
sampled version of tha.t of an analog filter. Let H(s) have 
simple poles and N, then 


C 

'(s) = = H(s) (s+dj|_)j 


s=~d . 


(2-2. t6) 


The analog f il-*.er iixu ulse response is given by 


h(t) - 2^ e ^ ^ 

. • i=1 

■ ^ -d,t 

h(nT) = 2^ ® ^ u__^(nT) 

i=1 


( 2 . 2 . 17 ) 


( 2 . 2 . 18 ) 


where T is the sampling period. 

N .1 c. 

h(nT) ^ ^d7T 

i=1 i^1 1-e z 


( 2 . 2 . 19 ) 


-di'r 

Obviouely, this uses the mapping, s+d^— '>■ 1-e a 

simple pole at -di. Since H(z~e‘’ ) is the Fourier Transform 

of the saiopled signal h(t). 

o<a 

H(e^“^^) = i i__ H(3(5-' 'F"wm) , where w = 2 t!/T. 

^ k=~ao ® ® 
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The equivalent analog filter must be band limited to the range 
(-n/T, VT) ; as T becomes larger the effect of aliasing becomes 
negligible and the digital and analog frequency responses 
become compp.rable. 

Bilinear Transformation : Using this technique, analog 
filters are converted to equivalent digital filters using a 
transformation [3] that maps the s-plane into the z*plane. 

The bilinear transformation is given by 


2 (i^, 


T 


( 2 . 2 . 20 ) 


1 + z 

in the z-plane. The LH-s plane and the RH s-plane are mapped 
inside and outside the unit circle in the z-plane respectively* 
On the unit circle z = 




{ 1-e‘ 


.jXt-T- 


( 1 + e 




2 4. /Xl-T^ 


( 2 . 2 . 21 ) 


where and are the corresponding analog and digital 
frequencies, and their ‘relation is highly nonlinear. There- 
fore this transformation can be used only when the analog 
filter response is piecewise constant. Though realizable, . 
stable analog filters are mapped to realizable digital filters, 
the nonlinear relation between and -A- leads to a highly 
warped digital filter response. 
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Matched z Transformation : The poles and zeros in the 

8-plane are mapped onto the poles and zeros in the z-plane, 

— 1 — aT 

where the mapping is defined by s+a — 1-z e whore T is 

\ 

the sampling period. The matched z-transf ormation is not 
suitable mapping when the analog filter has zeros .at centre 
frequencies greater tnan half the sampling frequency# as the 
position of the zero will be greatly aliased in the z-plane. 


Least Squares Error Method : 

Ass\une the HR filter’s transfer function is of the form 


k 1+a, z”" + b, z”" 

H(a) = A n 3 

k=1 H-Cjj.z'” + djj.z'”'^ 

io. 

Let the .desired frequency response be [H (e ^) , 

I 

where are not necessarily equally spaced, 

error in in- ^uency is Q(©) 


( 2 . 2 . 22 ) 


1=1 # 2 #3 # . . . )(.M j 
The sqmred 


M 


Q(©) = Z 
i^1 


H(e 


165 - 

H„(e 


( 2 . 2 . 23 ) 


where 


e 


«..# (2. 2*2-1) 


Q(‘^)^in Ip found using methods of calculus on the computej; fij 
provides an exhaustive et’'dy of FIR and IIR filter design 
techniques. 



2.3 QUANTIZATION EFFECTS IN DIGITAL SYSTEMS - 


The analysis of digital systems is based on the theory of 
difference equations with constant coefficients. The word length 
of any digital system ia finite. Therefore quantization of the 
input to the system and discretization of the coefficients ere 
compulsory. Quantization leads to approximation of the signal 
value by one of an appropriate set of discrete levels not from 
a continuous range. This nonlinearization makes an analysis 
of the system difficult except for simple cases. A statistical 
analysis may however be carried out to determine the resultant 
distortion. 

Quantization leads to rounding or ‘grouping' of input data- 
samples, the range of each sample is subdivided into 'interval! 
of equal width q, 

i-§-<x^ i+f 1=0, + 1, + 2... 

(2.3.1) 

The nonlinear quantization replaces each value of x by the 

j/'nearest interval centre *1-,'. The quantiser output 96^ can be 

* 

considered the sum of the input x and a round-off error 
referred to as 'quantisation noise* x^^ 

X* = X + X (2.3.2) 

The noise x may increase the variance of and bias the mean 
v' of the estimate x^ but for an important class of random process 
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it has been shown ttiat [it] the biasing effects of quantissa-, 
tion avei-age out approximately or can be predicted and 
corrected even with coarse quantization. The statistical 
theory of Amplitude Quantization was put forward by Bsnnet fiS] 
but is essentially due to Videow [ll]. His principle results 
f*re as follows : 

Let x,y be two random variables with probability density 
functions (p-d.f) p(x), p(y) and joint density function 

p(x,y). Let the characteristic functions be = E(e^®*) 

and = -ECexpC j (a^x +‘ a 2 y) ) ] . 

Tiu- quantizati.on theorems are defined as follows : 

1. If p(x) is 'bandlimibed ’ so that W„(a) vanishes for 

I * IX * r HI *1 

a j ^ - e (e > 0), then every existing mean E[x j is 

ccmpletel.v determined by E[(x')^] and the first ordejf probabi- 
lity den..s'ily function of x is xmiform between (-q/?, q/2). 

2. If p(x,y) is ' band limited ' so that ^ vanishes 

X ^ y 1 c 

for _ e, ^ - e, e > 0, then every 

existing moment E x^y^ (m,n = 1,2, ...) is given by 

l/(x')°' (y’ )^) and x and y are uniformly distributed end 

\ j tL 

sIb tisticalXy independent. ■ 

»row also gives a rigo^’ciig justification of the 
relations known as 'Sheppard's corrections'. 

B(x) =E(yM# B(x^) =B((x‘)^) - 


(2.3.3) 
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Watts [l4] has exteaied Widrow's theory to the- case of 
zation by shifted class intervals, i.e. 

% 

a, + i - q^/2 < x + i + q^/2 i 

t (2.3.4) 


+ 


k q. 


qV2 


< y a2 ^ ^2 '*■ 


The input-output characteristic of the quantizer (x* vs x) and 
its composition into shifters, gains and unit quantiser, are 
shown in Fig. (2.3.1). The p.d.f. p(x') of the output of the 
general quantiser is 

OtO 


P(x’ ) 


i 

r 


Z [ 

nr— iTsn 


P(x)fi [|^ 


c 


(| ” a)]dx].6(~ «C-Ii) 


(2.3.5) 

f^(x) is a window function s.t. f(x) = 1 for - x^ i, 
f(x) = 0 elsewhere. 6(x) is the Dirac delta function. 


The characteristic function (c.f) of x* is 


¥ 


. oo 

M = , Z e^P(-Sc<(o-a)) exp(-j2,.to) 


k=— oo 


i(ar~2Tik} 

(2.3.6) 

If X and y, with joint p.d.f p(x,y) are each quantised by" 

•» 

general quantisers with parameters (a^,q^,c^,r^ ) » i « 1,2, ..., 

the resulting c.f. of x' and y' is W , ,,(a,,a_), 

^ X ,y 1 2 

w,, ..(«■,.-,) = Z Z- exp( -3r.«, (c.-aj).exp(-ja2r2(c2*-a2)‘ 

^ $0 * ^ lrs=*^/jo X=«* ^ I I I I & c c. 


exp(-j2itka^) exp(.j 2 Tila 2 ) - |^, - |^)- 

sini (a^r^-2TCk) sin i(a 2 r 2 - 2 xl) 

?(a^r j-2iik) ’ t(a^2~2wl) 


(2.3.t) 



H 





Thia reei'lt may be extended * . my inamber of related , 

variables. Quantization has previously been described as e 

% 

addition of independent noise to the signal input to the 


quantiser, i.e 


X ’ = X f X 

q 


then E:(x':. ’) = E(x+x -y+y ) = E(xy + x y + xy + *^7-) (2.3**^) 


In t^-rins of the c.f. 
B((x’)'“(y')'^) = 


(-3) 


m+n 




1 


'^oc” ^0^ 


a 1=02=0 


(2.3.9) 


let the irrespective meo.ns be subtract^'d from x and y, and 
W = V/ (--i-ll - - 1^) be the c.f of p(x»x, y-y ) , V 

^ ^ 4. ^ 4 1 42 42 

is shown [l4, p 21 lj that if the c.f. ¥^(a^,a^) be zero for 


a 


i and 


n_ 1 ^ then 
2l qj 


E(-n ’y' ) = l/{xy) (2.3.10 

i.e,, tlif quantisation noises are uncorrelated with I *.h tb * 
input signals and with each other. 

Use of M tiler ; 


If a tnndom variable n oitisiies the first & sec'’o ^ 
quantizing theorems, and if x an^ >■ statistically i*-’' -. 

then the sum x+n also satisfies the quantizing the'^teins, 
less of the distribution of x. Chang and Moore [}6] stu’ttM fVy 
modified digital correlator shown Fig. 2,3 2. 



2,H 




?, ’ 1 


Let tije O', of + ^ nd n.(t) 5'e W («*>,) - W (a«) ree~ 

1 d. 

pective.ly. r^nng and Moore [l6j show tiiat t^*e ccrreloMon 
function of outputs of the quantizers y^'(t) and y 2 '(t) Pod 

ttot of the arbitrary inputs x^(t) and X 2 (t) are equal .Iff the 
auxiliary noise vnvr-forrcs n^(t), n 2 (t) are such that 


¥ (- 2-“-!:;) 
“i ‘>1 


k / 0 


and 


W I 

XI,-, 


2 Til 


) = 0 1 3 ^ 0 


(2.5.11 ) 


where q,,qj are the step size^ of the qtxantizers, i.e. 

] i. 


B(x^(t^) X 2 (t 2 )) = B(y|(t^) y2(’<=2^^ 


(2.5.12) 


Th' se conditions (eqn. (2.3.11)) imply that n^(t), n 2 (t) should 
be such that (i) they are zero mean (ii) they are independent 
of the nois'-' signals and of each other (iii) their p*d .fs 
shouJ ' be such that their characteristic functions have periodio 
zeros with period 27i/q^» Su/q^ no zero lying at the origin. 

It is shown [l6j that the rectangular distribution 


P(n) = 1/q |nj.< q/2 

- 0 elsewhere 

aaliaties eqns . (2.3.11). The p.d.f's which are m-fold convolu- 
tions ra = 1,2, ... of the rectangular distribution with itself 
also satisfy equation (2.3.11). 
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' The general form of the characteristic function of thtsj 
family of r.v.s. is 

oa ' _ m. 

Wj^(a) = [1 - i^rl ^ • F(«) (2.3.14) 

where are nonzero positive integers, and F(a) is s.t. 

F(a) 1, F{0) = 1. For the corresponding auxiliary noise 
n(t) to he zero mean, P’(0) = 0. 

% 

DEVICES FOR CORRELATION BSTIMTI.ON : Correlation techniques are 
of practical importance in technology. An important use of 
these techniques is the measurement of target range and velocity 
hy radar and sonar systems, and determination of the transfer 
functions of control systems. Correlation devices may be either 
analog or digital. The various categories are direct correlator, 
^digital correlator [l4j, stielfcoe's correlator [l4], modified 
digital (polarity coincidence) correlator [l7j» relay (modified 
stieltjea) correlator [tS], modified relay correlator. Fig. 
(2.3.3) shows self explamtory schematics of these correlation 
devices. 

The Signum Function and its Approximations 

The signtua function is a step function defined as 
follows : 


agn(x) 1 


X < 

0 


i 

1 “ ° 

X = 

0 

(2.3.15 ) 
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Thh uni' functiun i?. ■ ■fine'^ tiias 


!1 ( -'i. 't ,' = 0 X < D 

'j = J X =0 

= t X ' ’0 

mJ 

]n tern-- of '’'■'x), ■*‘l\eref o , 3^"{y) is "^iv 'u by 

ogn(x) = 2u(:.)-1 (''.3.1') 

'The unli ."tep u(y'' may be repres' the limit of ••■^ny 

uoritinumj , 


1 . U ( X ) 

= 1 irn <l 

u - ' 2 

+ ™ tan”^ (ax)] 

(3.d.l7) 

? . U ( X ) 

1 3 ra r 

a » oo 2 ‘- 

e,’f(ar) + s 

(b.3.l.'0 
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Tiiese fuiicti" ';s are illustrated xi FI j;. (^.3.i). 

?,i A. titat^stlcal Ooaiue-SarH. Led- J> ta Filtei'. 

In a recent pub3.iGa t i c»' f? ] a new technirue i i.. r the 
hardware implementation of digitvj filters i‘.' intrcducel. 
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(li) 


(iJi} 


(,lr) 



’{*)• "IC t«a"' I «? ) ? 3 



fig* 5^-.? tin Appryi-iaa lions of t.*'-*. Function 
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By of this teohni.'uie, tho con'j'lexi of the h. • '’..'ar© 
requir'.;d is rought to be reduce'^ The |;roiOse'’ fjitor 
requires only one summe-.', one shift register of Ifr'nrth cnuel 
to the order of the filter, either a simple RG low peso 
filter or a digital counter, and a generator of noise 
.vniformly distributed in amplitude and independent between 
samples. The principle of operation is based on polarity 
■ ''oidence correlation [S], Fig. (2.'3.3d), and is similar 
to dithering employed to over-come the effects of quantisation 
in video data for e.g. [9j,[l0]. It is used when lo\ input 
sample rate are sufficient as this filter requires a higher 
sampling rates than the direct. 

The proposed filter implementation technique is 
developed as follows;- 

Let be a random va,riable with probability density 

f nr oH An p(Xq) . 

ju ' 0 n(,t) be a noise process, stationary, whi^o, uniformly 
distributed in amplitude and independent between sctmples, 
lying in tho range -L'<n<C ¥.• 

The range of is included an that of n(t), or at moat 
the two lunges are equal. Assume teat the range of is 
equal to the range of n(t) and let it be normalised so that 

- 1 < ^ 1 


(2.4.1, ) 
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?hus , Ih'. T density of n is 


Pj^(«) =■- 1 ~1 <a < 1 

= 0 false;, 'har* 




HOi-/ consider the rani'.oi ''-'r- r 1?“ M c- z 'onoretr;; thus 


z = s, 2 n ( Xq f n ( t ) ) 
The conditional deneity p( z 


f > ■=; > 


x^) is p:iven by 


p( z|x^)=6(z~i ) .p(n( + 'S> + 1 )p(n(t)< ~x^) 


1 

6(s-.i) i j dn + r.('^'4-i) 

1 

J 


0 


in 


(2.4.4) 


I-!---':, 1-.., 

6 (Z-' ; + 'M ' + 1 ) — - 


Thfa' '‘cmii itioiV) 1 i'- c is C i-cT’ ir ,';4v ’O by 

E(z|x^) --- j 0 :/(-•( I V 


i f : 


J 

z 


I - ' (uM ) 


(2.4.4) 


( '.4 .7; 


I . hi 


I *T 


(?./■. ■''0 


r 2 ./. .'■}) 


S(Xo) 


S(i') ^ 

'J 

An intuitive evplon' tier of eon. {'•.'■.IC) r.o 
^'iven. Consider the case when x^ varies slowly in ti. 


f 


0 ) 
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w.r.t. n(t). The dithering aignel nCt) ie added to my 
N tlTies in a period of yiseccij each time the atim p ie formed. 


p = + n(t) 

and z = 8gn (p) 


We can see that if he deterministic and x^=0, because 
the probaV)ilities of n(t) bein^ '.ess than gre?’ than 
zero, are equal, therefore p shall be (+ve) or (-ve) with 
equal probability. This imp3i<=s that s'- is +1 or -1 with 
equal probability or its mean shall indeed be zero which is 
the value of in this case. If' x^ be -t-ve (-ve), tlie 
tendency of z will be to be +1 (-1 ) a larger number of times 
than ~1(+1). Avora.ging z over a large number of n(t) saraplGs, 
will yield some estTmate of 

A digital filter may be given by the difference equation 
N M 

-yk= Z (2.4.12) 

i=o j=1 

where is the input and y^^ the output. Consider the 
formulation 

IT ii 

= sgn [ 21 + 21 (2.4.1U 

i-0 3=1 

I 

where n is a noise process with density as in (2.4.2), 
stationary , white, and independent between samples, '’’hen 
by equation (2.4.10) 
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H 


S{y^)= 2_ ""A-i + 21 

3=1 


1=0 


(2.4.14) 


whore i=! assumed deterministic. Equotions (2.4.13) and 

(2.4.14) suggest a filter design, the output of x^hich after 
averaging many samples yields an estimate of the filtered 
u^. This is the quantised digital filter obtained from the 
linear or unquantised difference eqiiation given by eqn. 
(2.4.12). The output y^^ is followed by an averager, either 
an RC averager or a digital counter. The quantised filter 
obtained from the follovring difference equation 


k 


"k + '’gk-l "2 yk-2 


(2.4.15) 


is illustrated in Pig. (2.4.1). 

Therefore the output of the nuantized digital filter is 

an estimate of the true digital filter output. The output 

power spectral density is found b^^ ^-tranof orming tne discrete 

A OF R (m) of the output. 

«y 1/ 


H(z) 2 



] „ . z' 

3 3 



(2.4.16) 


is the ratio of the output and input spectral densities 
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Finite Averaging Birrors ; Since infinite time is not 
available for integration, approximate methods such as the 
RC averager ( a digital counter) must be used. Let y be 


scaled and shifted such that O^yxl and B(y) 


0 


Since 


y is ergodio 


y = 1 for p^ fraction of the time and 


y=0 for 1 ~Pq fraction of the time. Since the dither is 


independent between samples, the auto-correlation of y 1 

Jr |<T 
|dij>T 


> 


+ (p^-PoXi- 4*- ) 


= 

^o 


(2.4.17) 


f being the sarrplo and hold interval. 


by 


The variance of the output of the averager is given 


‘Z _ /„ yZ\ X A RC 




(2.4. 18) 


/ 1\ T 

(P -P^) • 

\ X o 0 


T 


2 RG 2 He 




( 2 . 4 . 19 ) 


It shall be shown in Chapter 3 ', th" t the- output v-ri'^nco 
may, be decreased by a change in the ‘rtructurc of tho quantized 
filter, using a digital counter es tja averager. 



CHAP TEH 3 


P‘*;V!’LOrMhijT OF AiJ N-BIT PEEuBA-’’' 


OOAIloE-BAI-il L: 


.ii- 


A /: i/i'ER 


In this chapter +!ie statistics oi the output ' first 
ani secon^l order linear digital filter are evaluatec? (dec 5.1) 
Klrlin’s coarse “ s'-M pls'i-d a ta filter (Sec. 2.4) incorpora tin^?: 
a digital counter aver? .^er in pla. ce of the RC averager is 
discussed in Sec. 3.2. The coarse - sampled-data filter is 
modified to yield the II~bit feedback coarse- sarnpled~data 
filter, dj '-cussed in sec, 3.3, the output statistics '"•f vdiich 
can be made to approximate closely those of the output of th '■ 
appropriate linear digital filter. 

1 The Linear Digital Filter.' 

A discrete time j^stci!' is an algori h’.m that cor'a-cts 
the input sequence u(nT) ^nto another ;;equenco (the out(''i + 
y(nT), the two sequences being related h'/ the oporator 0 

y(nT) =0 [u(nT)] (3.1^1) 

A linear time invariant system is characterized by its 
impulse response h(n) [l,pp. 13-1^']. The respujise to u(n) 
is given by the convolution relation 


Qqi ___ 

y(nT) = h(mT)u(n.-mT) 


(3.1.2) 



. We shell consiler causal, stabl--' LTI Fyt-tei.'.s. he inrut 
an-l output sequences to a LTI liiscrete system may be relate-^ by 
a constant coefficient f'^th Drier dii'f erenc"' equation of 
the form 


M _ M __ 

y(nT)= £ aj^u(n~kT)f- y(n-kT) n 15 , 0 (3.1.'^) 

k=o k=1 

I I , characterize the r • stem and bj^^j 0- The transfer 

function of the system represet 'ed.by (3.1.3) is 


^ -k 

T.. 

-.k 

I - 

, . k=1 

The frequency response of the system, H(e^^ ) is given by 


H(e^^ )= ii(z) 


z 

k=o 


• 0 “V ^k 


25=e^*** 1' 


M 

2:\( 

k=1 


(3.1.5) 




where ^ is in ratiians. 

We shall consi(3er the specific cases of the first and 
second order difference equation i.e. when M=1 and M=2« The 
input to the system u(nT) may be either deterministic or 
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random. For each of the two caep-s the statistics of the 
output of the system shell be deternined. 

3.1.1 The First Or’der F-ilter: 

The first order discrete time system is given by 
the difference equation 


y(kT)=a^^u(kT)+a^u(k-lT) + b^ y(k-l‘T) 


(3.1 .6) 


Consider the case vrhen a^=o for simplicity, jb^j <1 for 
stability, (3.1.6) reduces to 


y(kT) = u(kT)+b^ y(k~1T) 


(3.1 .7) 


Let the input to the system be u(nT) U(rT), U(rT) is the 
unit’ step fmction; by use of the recursive reletionship 
of eqn. (3.1.7), it is easily shown that 


^kT 


a 


o Z ‘■T 

r=o 


(3.1 .B) 


;(a) The Deterministic Input ; ^f^hen ii(kT) is determinist-ic, 
y(kT) is deterministic and from eqn. (3.1.8) 


B(yj5.j) 


k 

L 


®^^krr-T ) 


(3.1.9) 


r=o 


k 


21 ^k~rT 


r=o 


(3.1 .10) 
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A time autocorrelation function R (m) may be defined for 

' 

the -leterministic output .given by 

Ryim) = ISh'T JL ^kT-^k4¥T' (3*Uil) 

k=o 


Jim 

X.-S’jc 


K 

L 

k=o 


(a 


o 


r-'o 


Ut 


■rT^ 


k+m 

v 

Z- 

i=o 


U' 


1 k-}-m-iT 


) 


(7.1.1?) 


lim 


1 

K+1 


K 

k=o 


a 


k 


fi, 

'V»=r/' 


•"+"> ,rH, , 

.i- “l '^fT“k nn-iT J 
1=0 


(7.1.17) 


(b) The random input ; Let u^rn zero mean, white, stationary 
independent noise, specifically:- 


B(u^g-,) = ny^=o all n 

p(u^j,u— ^) = p(^jj:{:57p) p(u^,^, 


Than the expected value of the output using the oqn. 

(7.1.8) 

^ r 

E(yi.T) = J. l>t “WT 3 

= •£. ’'i' 

** 3^=0 


= 0 


(7.1 .1»'l 
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The ensemble autocorrelation of y{kT) is 


R (tok+m) = B(y^j yp— j) 


» k k+m 

= "o I., Z ^ 

r=o 1=0 


r+i 


Efu 


U: 


k-rT ‘^k+a-iT 


_ _ k k+m 

R (k, k+m) = T Z <5(ni+r~i)b . 

y u u. 1=0 


r+i 


1 

(3.1.16) 

(3.1.17) 


The autocorrelation of may also be derived as follows; 


E(y(kT)y(Ic+mT))=a^ E(u(kT)u(F+mT) )+aQb ^ rs(u(kT)y (Jc+mTfT) 

+i!:(y^1T ) . uCk+lnT)) ]+b^ E ( y (F^T t ) y ( Km-TT ) ) (3.1.18) 

The cross correlation terms P (k^l, k+m), R (k,k+m~1 ) 

y ,u u-»o' 

nre first evaluated 


(i) ■'‘y.u'*"-'’ "k+m^ 

m > -1 R,^ (k-l,k+m) = 0 (‘^.1.19) 

This is because y ( F~l'T ) is a function only of u(F--lT) and 
all the previous values of u, and is independent of the 
future values of the input. 


m 




= ®[ (®o“k-l'''’’l ^k-2^ “k-1^ 

j2 


a 0 ~ + 0 
o u 


(^.1 . 20 ) 
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n» <; -1 




= 0 . a „ a ^ b , l"-'l 


(, 


(ii) R^^y(k,k+m~l)= E[u,^ 
m < 1 


![uvy 


k-f m - 1 


] = 0 


m = 1 


.2 


= a „ a 
o u 


( 




m > t 


B[^i 7 ) , 1 ] = a 

k/ k+m- 1 o u 1 


Therefore , 

V 


(" 


R (k,k+m) = a^ (l+h^+(b^)^ + ... + m--R) 

y Oil 1 ^ 


= a^ b“ [ (H-b^+b^+ ... + b^"^).o'^] m>0 

- =a2 a2 [ u - h ^^; ... b," 

^2 , >,2^+2 

B , 1 — o - P 

R (k,k+u,' '• . tyf m = 0 

y 4 ^ 


, 2 ,^ 4 , 


’■”^1 

, _, 2 k + 2 ' 

.2 „2 ..m • 

j cr b. ' ^ " "' 'A^ 

o U 1 . 2 


>0 


l-b 




■.1 ,21 ) 

3 . 1 . 22 ) 

!.1 .23) 

i.l .24) 

ra < 0 

(■^.1.25) 
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T'l* ini'Ut to tha nysi^m is U(jT), U(jT) bein/? the 
unit otap function. The process is necessartly nonststion 
r gnrlleuK of the nature of Hence the outrut of th^* 

nyster y^^ is nonstationhry , which fact is ratified 

by eqn. {' 5 .]. 2b). '''fiit-n the systc-T! reaches steady state j.' 


Ry(Ki) 


o2 ^2 

B O 
0 U 

1 h2 ' 

1-b, 


» t 


with k very large, and |h^| < 1 » eqn. (5.1.25) tends to 

m = 0 

I 

(5.1 .26 j 
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O il 


m 


1-b 


1 Imj 1 


In steady state the power spectrupi of th'- output is 




2 ? 2 2 

a a 

0 u ^ 

1 - b^ 1-b^ 


• j '*> m 


(3.1 .27) 


21 c-J'-"" i,'!' 

i71=^1 


) = 


2 2 
af a 


° ^ [ 1 12 Z coat-nij ] 


( 5.1 


1- b 


1 


m=1 


Th value for S„(‘*^) ob taint'd usin;, the reJntion 5y('^) 


y 


u 


I(z) n(z'"'') _ fora two side-' input - n'! th't ni to in 'n’ 

2i”3 


in eqn. (3.1^6) are equal. In "tf- ;y 
expression for S (<^^) mei' be used. 


•-;l -1. 


uiO ucup3. 
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3.1.^ Th.' Second Order Filterj 

The second ordor linear system is defined by the 
diff.-'rence eqiiation 


y(kl') = u(kT)+ b^ y(k-lT) + b^ y(’--2T) (3.1.P9) 


The transfer function of the system is gi^'en by 

a 

H(z) = 


1 “b ^ z 


•1 


(3.1 .30) 


bpZ 


The inverse-z transform of Fl(z) gives the impulse 
response h(nT) of the system. y(kT) is given by the 
convolution integral 


<50 

V 


yCkT) = l___ u(mT) h(lc^T) 

■ ’ m=o 

H(z) is given by eqn. (3.1.30) 


(3-1.31) 


H(z) 

= ^0 

d-.p^z’*'' . 

(I-P 2 z“^ 



b 

where = — 

2 

+4bjV''^ b 

, P 2 “ 

-(b2+4bp) 
^~'2 ^ 

1/2 

(3.1 .:»?) 
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( 1+P^ 

~1 , „2 ~2. 3 - 

Z +P^Z +p^z 

.3 

+ ♦••)( 1 f z 

■ 1 ? -.? 3 --3 \ 

+P 2 Z '+P 22 +...) 
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i 

1 — 1 

P2~^ 

]) 

(5.1.35) 



n=1 1=0 




Hence h(nT) 

= 0 

n < 0 



(3.1.34) 


= a 
^0 

^ ^i n-i 

,ZL— Pi Po 
i=o 

n ^ 0 
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T he Dotbrm in iBt io '(k/ pv3 of eqo, (3»i.3l) 

ives boUi the output of the Byrf ira an.'^ its expect'''^ 
value, wh 'n u(^'’'') is d etermini stic.The time autocorrelation 


is <Jlven by R„(m) 

A 


H,. (in) 


1 ' 'n 1 

K-^oo E+T 


i y(hT) y(ra’) 
k=o 


lim 

K '^ao 

gO cjC 


1 


L T. L h(lT)u(k~tT)h( jT)u(Ff'^T) 
k=o i=o 0=0 




iPO 


= L L h(iT)h(oT) JiilU'KTT A «(k^T)u(trr'TT) 

i=o 0=0 . k=o 

Now if the viine an tocorre.lr.tion function of u(kT) is 

R (m) = ^I'r- Z. u(kT) ’-(kFnT) 

u' ' K+1 

substituting this relation in ('^»1.35) we obtain 


"O cK> 

R (m) = Z. L h(iT)h(oT) .R^(iii~ 0 +i) 

y' i=o o^"*^ 


(3.U36) 


(^) T he Ra nd om Input ; to tn pyatcm is u(nT) U(oT), -ore 
u(nT) is defined by equ. (3.1.14) and U(oT) ia the unit step 
function. Prom eqn. (3.U3l) it is deduc ' th-. I 

S(y(kT) )= L[ L u(mT)h((E:mT) ] (3.1.37) 

rn=o 

eo 

-- Z E[u(mT,i h(k-mT) =0 (3.1.38) 

m=o 



The autocorifslri tion of y(kT) is defimod as 


R (k,k+m) B[ y(kT)y(k'}niT) ] 


<30 

r 


.50 


ZL li(T^T)h(E?i^T)0^ ^(n-J) 
n-o J=o 


«7<0 

.2 Y" 


u 


h(k~nT) h(k + n*.nT) 


(3.t.39) 

(3.I.4C) 

(3.1.41) 


n=o 


nary R (k,k+ro) 


Obtaining the exprension for the no; ■> > > 
by using eqns. (3.1.41) and (3.1, .34) setting k ■-»<«> for sterdy 
stat' onditionrj and obtaining the DPT of the etationary 


H (rp) is unnecessary as the relation 




H(z) . 

^ z=e'^ 




(3.1. • 1 


tC obtain the power spectrum for a two sided input yields 
1:he same result. Therefore 


S (.u) = 


2 2 
a at 
u o 


■y ' b e"='" ) 


.i2'>> 


1 


(3.1.43) 


2 2 
'^u * ^0 


14-b^+b2-”2(b^ t-b^b 2 -)cosCA'~ 2 b 2 cos 


(3.1.44) 

Tbs autocorrelation function is given by the inverse-z 


transform of Sy(z) 


W 8y(ml)=j^ j'sy(z) 


z d z 


(3.1.45) 



5*11 


for stable H(z), all its polts li'‘ within the unit e*rcle 


« = 


j2iifT 


f") -:3?j3tH on the unit- circle anl 


Ry{mT) is .::i- t, ty 




j-> 

2T 


■/ 


'•( . exp (jKfTk)* rif 


{3.t .?f’) 


1 

2T 


Tho eecont order difference equatiei'.i may bo simplified by 

% 

setting b^=0, in which case the .r'''cursive relationship is 


y (IT) 


u(kT) + b..y(k--21') 




The transf-='r func.tion of the systo’ii is 






1 


( 1-p^ 7,“ } ( 1~Pp7” ' ) 




;; f ± p 


ih. 


(3.1.48) 


The steady state power •’poctrum of the output is 


S (to) = -t. 

^ o' /z:rpl‘ ^ 


o^-E(r.) P(z"^)i 


2 2 
a a 
u o 


1+b2 ■" 2b2COs2<^ 


S (lo ) = 

y 


13.1.49) 
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Using eqn. (3.1.47), y{kT) may be written in terms of the 
input v??lu?a (u(iT), i=0, 1 , , . ' iff or k=2n, 

I 

as follows: 

y(kT)= a^u(kT)+a b^u(2^T)+ . . .+ a bJuCkrai^) 

U VC O c 

+ ... + ^ 2 even k (3.1.5C^) 

= aj^u(kT)+a^b2u(k^T) + . . .4-a^b2 u(F-r2pr) 

+ ,*, + bp^ u(1) for odri k (3.1«50b) 

Thus y(kT), k is odd (even), is a weighted sum of 
the odd (even) input values. 


Evalus-tion of the Autocorrelation function R (k,kim) of 
the output y(kT) requires prior knowledge of the cross- 
correlation “functions viz. R (k,k+m-2) and R (k-2,k+w;. 
Using eqn. (3.1.50) 


r 

\,y 

t 


= B[u(kT)y(k+m-2T) ] 


0 

2 / •» \ iii*“ t 

a^a (b^)~t 
o u' 2 2 


m < 2 ,m= all Odd (+ve) 
integers 

m= all even(+ve) integers. 



(3.1.51 ) 

(k-2,k+m)-BCy(F::2T)u(S+mT)) 

= 0 m>-2,m=odd (;-ve) integers 

~ m= all even(-ve) integers 


(3.1.52) 
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Now 

]fty(k,k+m)=a2 E[u(kT)u(r+mT)]+a^b2 B[y (i::5T)u(S+mT) J 

+ ^0^2 ^ f ^ ^ T ) y ( F+'m-'2 T ) B [y ( k-2 T ) y ( &m-2 T ) ] 

(3.1.5?) 

Substituting (3.1.51) and (3.1.52) in the recursive eqn. 
(3.1.53), for k=2n+1 or k=2n and k+rit=2r or k+m=2r+1 


Ry(k,k+m)= cr^ ( 1 +b 2 +'b 2 + . . .+ (h 2 )^) . in=0 

2n+2 


-2.^4, 


.2 ^n^ 


2 _2 J "’^2 


(3.1.54a) 


= a 


Qb 2 .a^ 0 ^b 2 “^^"’^[l+b 5 +h^+...+ (b^)^ ] rn=?,4,6,... 


= ^2 af , b ”''2 1 :^% 


2 2 
2n+2 


(3.1.54b) 




1 v2r+2 

2 2 , '{ m {/2 r ^”'^2 

‘o‘’u-’'2' t- 


1-b; 


0 I ni| — 1,2,5, ... 


m=-2,-4,-"6 , . . 

(3.1.54c) 

(3.1.54d) 


In steady state conditions when k Is large, the Autocorrelation 
is stationary i.e. 
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Ry(rn) 


1 im 
k— 


Ry(k,k4'm) 


2 2 
a 

o u 


1-b 


2 2 
a a 
o u 


1 -.b^ 


m = 0 

b H /2 
2 


i 


2 , 4 , 6,0 


1:^.1 


= 0 elsevrhere 

3.2 The Coarse-Sampled-iiata Filter. 

In (Chap, 2.4), the basic cob ta filtpr 

using an RC averager, was studied, ilerp we study the c''arse 

* 

sampled -da ta filter using a dlgita.1 counter- a/errger 1'O'f’ 

.< '.'’rrison with the linear digjtol 1'il‘o'' an' for the 
subsequent comparison of its output i tfi ilir oui.pii^ of itr; 
modification. 


3.2.1 The First Order Coarsc-sarnpled -iP. i i ! '• r 

The difference equation chare etc -"'i r in ''irfit 

order linear digital filter is 


y(lT) = u(lT) + b^y(l-fr) ('.2.1' 

The filter is replaced by the coarse- samj.lr'-i — •-* 1 " ' Itur 

which serves to reduce the comple:'-! i*:/ of thf hi^ircwaro 
realisation of the filter. The relation het-jcrn t'-: ;ln t 
u(lT) and the output y(lT) of the coco’se-sampled-dc ta filter 


is given by the equations 





z(Nl+k)= agn[a^u(lT)+n(Nl+k)+b^z(Nl+k- 1 )^ (3.2,2a) 

y(lT) = i 21 z{Nl+k) • (3.2.2b) 

" k=1 

A schematic diagram of the* filter is shown in Pig. (3.2.1). 
Referring to the figure, let u(lT), the input to the coarse » 
sampled-data filter, be defined as in eqn. (3.1.14). The 
input and output sampling periods of the coarse~sanipled - 
data filter are T secs. T should be chosen such that it 
is commensurate with the Nyquist rate of the input. The 
input should also be effectively constant over the sampling 
period T. The output of the ZOH is 


u(ia+k)=^ u(Nl+kT' )= u(lT) for 1«k^N (3.2.3) 

The dithering noise n(Nl+k) added is such that it is 
uniformly distributed in the range (-1,1), independent 
between samples, and white. It is normalized, so that 


C^o u(lT)] + |b^j 4 ^ 
and [a„ ^ 


(3.2.4) 


The output of the comparator o.r sign detector is 
+ 1 (-I) as its input, viz, a^u(lT)d-n(Nl+E )+b^ z(Nl+K-1 ) , is 
greater than (less than) zero. The dels.y of the feed 

back register is T’, which is the period of the dither 











n(Nl+k). The averager is a iigital counter that is reset 
'to zero every IT secs, 1=0,1,.,. It forms the weighted 

, t c 

sum -4— / z(Nl+k), up counting (down counting) by one 

1^1 

as z(Nl+K) is +1(~1), effectively averaging out the effect 
of the added dithering noise n(Nl+k). The initial value of 
the delay line, z(-N+N), (in this case it is e single 
element corresponding to T') is +1 with equal probability. 

gives the input-output relationship of 
the coarse -sampled -data filter. The statistics of the 
output y(lT) are evaltia.ted for 1=0, 1,2,... for the case when 
the input, u(lT) is deterministic and random. 

The Deterministic Input ; The range of the d eterrainistic 
input to the system u(nT) is such that eqn. (3.2.4) holds. 

The output' y( IT) is given by eqn, {3.^?.2). 

(i) The probability of z(Jfl-t-k) : The single element z(N1.4-k-1 ) 

f' 

is initialised at l=o and k=1 , so that the probability 
of the initial value z(o)=z(N-N') being +1 or -1 is 
half i.e* 

p . ( i )=4 = P , (-0 0 . 2 . 5 ) 

Prom eqn. (3.2.2a), for 0, k=1,2,...,N, 

z(Nl+k)=1 ^ a^uClT)+h(Il+k)+b^z(Nl+k-l) > 0 

n(Hl+k)> .e^u(lT)-b^z(IIH-k-l) (3.2.6) 


*• 
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The conditional probability I . (t ) 

^•^l+k { ^a+k-1 f 

la given by 

‘’wi+k I qa+k-1 >> -?„u(lT)-b,z(NH-k~0 


yp{n).dn 

-a ( IT ) ~b ^ 2 { Nl+1:- 1 ) 


1 -Fa^u ( IT )+b . 2 ( IIl+k~ t ) 


(3.2.7) 

( 3 . 2 . 8 ) 

(3.2.9) 


Therefore 
P 


J (-lj2(Nl+k~l)= 

^Nl+k ^Nl+k-l I 


1 -a^tt ( IT ) ^ 2 ( I'Tl+k- 1 ) 


(3.2.'10) 


Since 2 (lil+k-l) = + the four conditional densities 
may be evaliia.ted from eqns, (3.2.9), (3.2.10). 


fP 


z 


'ifl+k 


,(i)= I 


(1| 2(Nl+k-1 ))[p 


^Hl+k-1 ""i^i+kpNl+k-.l i 


2 


'Nl+k~1 


1+a^u(lT)+b^ 


(2(Nl+k..1 ))] 

1+a u(lT)-b, 

p (1) + -^-^ 1 

^Fl+k-l ^ 


(3.2.11) 


* ^,^^l+k~1 ^ ^ 


‘‘Hl+k 


( 1 ) 


1+a u(lT)-.b 


1 


+ bl P„ ( t ) 

^Hl+k-l 


(3.2.120) 
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From the recursive relationship of eqn, (3.2.12a) and from 
the initialisation z(-N4-N)= +1 with equal probability, 
the probability of z(ITl+k) may be determined. Stmll»>rly 


z 


Nl+k 


(- 1 ) 


u(lT)-b 


1 ’ 


+ b, p 
1 ^z. 


‘Nl+k-1 


(- 1 ) 


(3.2.12b) 


(ii) The Expected value of y(lT); 

y(lT) is a function of both the input u(lT) and the added 
dither n(Nl+k) as in eqn. (3.2.2). 

N 


y(iT) 


“ I Jr, 


A 1 N 

E{y(lT))=i J E[z(Nl+k)] (3.2.13) 

" k=1 

From eqn. • (3.2.2a), 

E( z(Kl+k))=E[sgn(a^u(lT)+n(Fl+k)+b^ z(I'Tl+k-1 ) ] (3.2.14) 

=- Z- • [ / [sgii(a u(lT)+n(Kl+k)+b. z(IIl+k-1 ) ) ]. 

» Nl+k-1 n ■ ' 

p^(n) dn] . p(z(Nl+k-l)) (3.2.15) 

Invoking the results shown in (Chap. 2.3), (3.2.14) reduces to 
B[z(Nl+k)]=E[a u(lT)+b. z(Nl+k-4),] 

U I / 

= 8 .^ u(lT)+ b^ E[z(Nl+k-1)] 

=a^ u(lT) [ 1+b^+b^+ . . .+ b^”^ ] + 


(3.2.16) 



K 

h. [a uCrrfT) ( ^ ' )~a^u{ l^T) .bf ( 1+b, f-b'.'i-, ^ , 0 ""'^ ) 

10 s O I I ] 1 


■i- ^ a^u(r-rT} 


±i{ ^ 'i ^ , 'U . 




•f a. 


u(oT) 


1-b 


®0 fib 


k 

L_ + a„ b^ 


{ 1 'h’”’^’ 1 

0 1 l-b^ * Jn... 1 


r=l 


(3.5.17) 


Substituting (3.2.17) in (3.2.13) we obtain 


N 


E(y(lT)) = I z F(z:Nl+k)) 

k=1 


a^u(lT) a^u(lT) N 

- - t^f;; 


(1-b^b 


J, ^1 ^ ^o TTTFp 


r j b^Ur-n u(Tirl’)>[i b^ ]J 
r=1 1 k=1 ' 


a^udT) b^a^ u(lT) (l~b^^) a^ 

1 —b j ,r ”/ * -u \ 2 N ‘ 


(i--bp' 


7 uCTTFr)] 

^■”‘^1 35=1 ' 


(3.2.18) 


For effective averaging, the length IJ of the counter must be 

with 

very large. Theref oreZjf-»-<« , from (3.2.18) we obtain 

N 


B(y(ll))= K?« W Z a(2(Ifl+i:)) 


1-b, 


k=1 


(3.2.19) 





By comparing eq. (5.2.19) and eq. (3.1.9), it is stt-n that 
the expected valup of tiie output of the coarse-sanjpled" data 
filter does not approximate the output of the filter it is 
required to replac** 

(ii) The Time Autocorrelation of y(lT): 


It is required to check the matching of the time 
autocorrelation Tunotions R^(m), RA(m) of the outputs of the 
linear digital and the coarse -sampled -data filters respectively 
let ft*(m) he defined as follows: 


p _ lim 

1 

I+T 

K 

Z 

k=o 

y(kT) 

y(k+mT) 

(3.2.20) 


1 

r+ 1 

K 

1 

■— iniwijpw 

V 

N 

I 

i=1 

N 

z ( Nk+i ) z ( lfic+m+ j ) 







(3.2.21 ) 

r 

_ lim 

1 

WT 

K 

L 

k=o 

1 

N 

i^ 

N 

V 

/ 

3=1 

Sgn(e^u(kT)+ 

n(Nk+i)+h^ 

z(Nk+i-;-1 ) ) .Sgn(a^ 

u(E¥T)-i-n(Nk+i+'j’) 


+h^z(Nk+m+3-.1 )) (5.2.22) 


Now since y(k) is a function of the added dither, the 
expected value of R*(m) is shown to he 


B 


(R^(m)) 


lim 1 X 

K-.o.i:Tr 




N 

^ B( z(Nk+i) .z(NE+m+ 3 )) 

3=1 

(5.2.25) 



When 15=0, for a finite N,E[R (a)] r^.'^uoes to 


N N T 

+ 7 ) S( z(Nk+i) 2(Nk:+ J ) J 

.0^ 

lim 1 5 N N 

K^., CN-^ -I, _LB(z(Nk-.i) 


I 1 


N 


7 


zdIK+j))] 


(3.2.24) 


Using eqn. (3.2.17) and the results derived in (Chap. 2.3), 
,(3^2.24) reduces to the recursive equation; 

E{R^(o)) = L ^ [U"-N)a2 uRkl')+2a^b (N-.1) 


NaQu(kT) 


a 


1 , 

u(kT). ^ +ao( — -f. L u(E::?T)’j 

° l-bt ® l^h, r=1 ^ '> 


+ b 


? i 

« jf — 4 


I 

L s[z( 
i=1 0=1 

In the limiting case, as N->toO , 


^(Nk+0^1)]] (3,2.25) 

,Rj) tends to 


Mi“ ^ J, 1, ‘ -i-' > ] ] 


N U 


(3.2.26) 





Comparing Eqn. (3.2.26) with Sqn. C^.I.O) fi(RA(k,k)) 

Hi — O i 1 c»0 y 

as given by (3.2.26) is see'i to be greater than RA(k,k) 

y 

evaluated for the linear digital filter. Similarly 
forjmf^, from eqn. (3.2.23) and using (3.2.17) 


K 




1 


k=o r 


N N 
i^1 j=1 


(z(Nk+i)z(NF+m+j)) 


dr ^ "V rN2a2u(kT)u(lj:+mT)+N a b u(l+m'T).[ f. 
K+] r ° ® ' i=1 


B(z(IIk+i-1 ))}fN a b.u(kT) T. B(z(Nl?i-m+ j-1 ) ) 
J o 1 


N jsr 


+ b‘ 


•• '7 Z B(z(Hk+i-.1 )z(Nk+mt- 3-1 ) ) ] 
! j =1 0=1 


(3.2.27) 


dr a qU ( k'P ) u ( Ic+mT ) tK a^b ^u(lc+mT), 

p(z(Nk-.1+N))+ Z E(z(Nk+i-1 ) )| + llap^u(kT). 

N - j 2 I 

jE(z(Fk+m~1+N))+ Z B( z(Fi:i-in+o-l! >>] + I Z 


E(z(Nk+i-1 )z(Brk+ra+o-1 ) ) ] 


i=1 r=1 

(3.2.28) 
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K 


*1 iX 

ttxT ^ C ( kl ) u ( k+mT) +Na _ b . u ( kSnT )• 

ixi j Q ■*vj^ ^ Of 


qaoUfE^T). ^ . |-’ ,tl(r-, 

+ I (N-1 )a^ 


l-b 


N 


k-1 


U(k0?) ^t~^i 


l-b 


1-b 


+ a 


1 , , N 

i-b, bj-b, 


o 1 ~b 


1-b 

-b 


k 1 

^ tN(r-l) kTFT )j j+Ha^b ^ u ( kT ) [ ^^ a^u ( lc+m-1 T ) . 


+ a. 


IT 

l-b. k+m-1 -v , a^u( 

'” ud+m-T-r'Dj +j (N-1 ) -2.p 


'1 

a u(.k+mT) 


a u(k+mT) b.-bf 

-2-5 i . 4^ + a„ 

1 -b^ o 


1-b 


1 


'1 

N N 


l-b^ b,-b® 

^ ~ 4, ' 


u(k+m-rT) j ]+b^ £ E(z(Nk+i-1 ) z(l!Tk+m+ j-1 ) ) ] 

(3.2.29) 

In the limit, as N~ 3 *»o » B tends to 


K 


2aib, 

o* 1 


dr )utk+mT)+ . 


u(kT)u(k+mT)+b^ T ^ X B(z(Nk+i’-l) 

N i=1 .1=1 


N H 


i=1 1=1 


z(lTk+m+ 3 - 1 ) ) ] 


(3.2.30) 





Bqn. (3.2.30) is a recursive equation; comparing it with 
(3.1.13) for m?^o , we see it yields extraneous terms 
apart from those encountered in eqn. (3.1.13). 

(b) The Random Input: 

The random input u(lT) to the system, with density 
function p^(u), is defined by eqn. (3.1.14) and its range 
satisfies eqn. (3.2.4). Eqn. (3.2.2) the input ~ 

output relationship of the filter; The statistics of the 
output may be evaluated as in the case when the input 
is deterministic. 

(i) The conditional probability of z(Nltk) : z(lTltk), given by 

eqn.’ (3.2.2a) is a function of both z(Nl+k~l) and u(lT). 

The probability p„ (l) and p^. (-1) are both equal to half. 

^0 o 

• Therefore the probability of z(Nl+k), conditioned on u(lT) 

and z(Nltk-l), may be evaluated. 

z(Nl+K) == Sgn (u (1T)~- + b^z(Nl+Ei^1)+ n(Hl+k)) 

Therefore, 


P 


^Rl+k 


^ (lj z’(Nltk-1),u(lT)) 

Nl+h-i ,u(lT) 1 


1+a^u(lT)+b^z(Ml+k-1 ) 


2 


(3.2.3la) 
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z 


Nl+k 


^Nl+k-l ,u(lT) 


(-1 


z(Nl+k-1 ) ,u(lT) ) 


1 -a^u(lT)--b^ z(Nl-<-k-1 ) 


; ■.2.3tb) 


(ii) The Expected value of y(lT): 

E(y(lT)) = i 21 E(z(Bl4-k)) 

k=l 

Expected value of z(Nl+k) is given, by 


E(Z(Nl+k) )=.S[sgn(a^u(lT)+n(ITl+k)+b^ z(Nl+k-1 ) ] 
=E[aQu(lT)+b^ z(Nl+k-.l)] 
=aQB(u{lT))+b^ B(z(Nl+k-.1 )) 


(3.2.3' 


By analogy w.ith eqn. (3.2.17), (3.2.32 ) reduces to 

1, (1-bf) 1 ^ 

B(z(El+k)>= a^ B(u(lT))+a^b; ^ ^ 

1 1 r=1 


^N(r-I)g(u(i^^,j))j 


(3. 2. '53) 


Therefore , 


. N 1 -b; 

1 V 1 


1 

k '-'=1 


B(y(ll))= ^ [a„ 




U 


1 

■c 

r=1 


t,N(r- 1 ) , 

•'Oi 


^o^u ®'o^u 1 


1-b. 


1~b^ * F 


. b 




1 1-b 


1 


F 

F 


^■*^1 / V F(r-1)^ , -j 

• V =1 ^ ^ ^ ^ 


(3.2.34) 
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In the limit, as E(y(lT)) tends to 

B(y(lT)) = ^ ■ (’.2.55) 

Therefore, when N tends to infinity, E(y{lT)) agrees with 
■B(y(lT)), (eqn. 3.1.l4h) if T)^=0. 

(iii) The Ensemble Autocorrelation function of y(lT): 

The autocorrelation ftmction of y{lT) Is given by 

H^(k,kfm)=B[y(kT)y(l+mT)] • (3.2.36) 

. N N 

= •^7 J E(z (Fk+i)z(lTk+m+j)) 

N i^^1 j=T 

Using eqn. (3.2.2a) 

. N IT 

M(k,k+m )=' Y Y. B[sgn(a u(kT)+n(lTk+i)+b.Z(lTk+i~1 )). 

^ i=1 3 = 1 ° ^ 


sgn(aQ u(k+mT)+n(ITk+m+ 3 )+b^ (uSi+m+j-l ) ) ] (3.2.37) 

t 

For m=oi“ 


IT ‘ N IT 

R* (k,k)‘= B(z(Fk+i) + Z. E[sgn(a u(kT)+n(Mk+i) 

y i=1 i=1 3=1 


+l)^Z(ITk+i*.l ) ) , sgn(aQu(kT)+n(ITk+3 )+b^ Z(Nk+3~1 ) ) ] 
, IT IT 

[K+ Y. "T B[ (a_u(kT)4-b., z(ITk+i-1 ) ) (a u(kT) 
i=1 3=1 

+b^z(rTk+3-1 ))] 


(3.2.38) 





op 
0 “ ^ 


N 


K fN+(N^-N)a^ b f E(u(kT)z(I[k+j-1 ) )4-(:;-.i) . 

N \ U C’ 1 


M 9 ±L - 

Jr^ E(u(kT)z(lTk+i-.1 ))+b^ )_ 2- 3( 2 (iTk+i- 1 ) z(im+o-1 ) ) j 


iT r 


i=1 3=1 
i?^3 


(3.2.39) 

It is required to evaluate the cross -correlation function 


(kT,KE+I...i'r')=E[u(kT)z(NkH- 3 -l)] ,j = 1,2...,N (3.2.40) 

for j=l 

B[u(kT)z(Nk+ 3 ~ 1 ) ]=E[u(kT) z(k%-T+N) ] (3.2.41) 


From the structure of the filter, u(kT) and z ( NE^f+N ) 
are independent^ the, latter being a function of u(k^T), 
u(k-2T) ,,.u(oT) , of which u(kT) is independent. By definitiot* 
B(u(kT))=0 hence (3.2.41) re4uces to 

r 

^ E[u(kT)Z(N¥rr+N)] = B(u(kT)).B[Z.(5rlcrf+N))]r=o (3.2.42) 
for 2 4 34 IT 


\ 4 -f 

B[u(kT)z(Nk+ 3 - 1 ) ]=B[u(kT)sgn‘(a^u(kT)+n(Nk+ j-1 )+b^ z(Nk+3-.2)Ji'’] 

=E[u(kT) (aQu(kT)+b^ z(lIk+3-2) ] 

=ao o^+b^ E[u(kT).z(lTk+ 3 - 2 )] 

^o ^u [ 1+1^1+^^+* ] 


2 4 34 N 


(3.2.43) 



Substituting eqn. (5.2.43) in em . (3.2.39) we obtain 


H-;.(k,k) = [N+(II^wN} .a^ cr^ 2 (N-1 )a„b , £. a.o^. 


1~b 


j-l 


N 


o u 


^ w ^ CA. _ VJ- • TT 

o ] ® ^1 


P N 

+b^ r ^ l[a(lTk+i~1 )a(llK^ j-i ) ]] 

i^ j = 1 
i?^3 


(3.2.44) 


- JL 

"r 


[H+ (K^-]?)af 0^2+2 (T ,.Oa' 




o u 


O 1 u 1 -i 


b . -b 


t ^ (l-b,)' 


+ b 


P N K 

1 "L L E[z(llk+i~1 )z(i »3-r’M 
i^ 3=1 
i?^3 


(3.2.45) 


In the limiting case, when F-^, R. ,k) tends to 

y 


lim 


N N 


.L .1 -.U(Nk+i-l)z(Nk4-3-.l)] +' 

IT IL 1 ^ 1 

2a^0^b^ (l-b^)“^ (3.2.46) 

Comparing eqn. (5.2. '6) with eqn. (3.1 - W) for m=o, we see that 
(3.2.46) contains other terms, apart fi ‘fu those contained 
in eqn. (3.1.17). 


For m^o:- 


1 ^ 

R^(k,k+m)= J-sT T J B[sgn(a u(kT) I n(Nk+i)+b, z(Nk+i-1 ) . 
^ i=1 3^ ° ^ 


sgn(aQu(k+mT) » »>(Fk+m+3 )+b^ z(Ih; I'H-j-l ) ) ] 


(3.2.47) 



Therefore n{Nk+i) and n(Nk+m+3) are indep' 'ident hence 


N 11 


l^(k,k+in)= 2^ Z~ BCaQU.(kT)u(k-rinT)+a^b ^u(k+mT) z(K^k+i~1 ) 

' N i— 1 n~1 


+a^b^u(kT) z(Nk4-m+j-1 )+b^ z(Hk+i-1 ) z(Nk+in+3-1 ) ] ’ (j.2.48) 


N N 

w [a b, .N. T B(u(l+mT) z(lTk+i-1 ) )+I?a b T B(u(kT). 
2 o 1 o 1 




z(lTk+m+j-l)+b, T T E(z(Nk4-l-.l)z(Nk+m+j-1 ))] 

^ i=1 3=1 

(3.2.49) 

To evaluate Ra (k,k+in), therefore, the cross-correlation 

v - 

functions E(u(”k4-mT) z(Fk+i-1 ) , ■E(u(kT)z(Nic+m+3-1 ) ,i=1 , . . . ,N, 
3 = 1 ,..., N need to be eva.lua.ted. 

B [.u( Ic+mT ) z ( Fk+i- 1 ) ] , i=1 , 2 , . . . , N; mj4o 


''B[u(k+fiiT) z(Fk+i-l)]= 0 m > 1 , for all i 


a bi-'.cV'U-/-') . 

o 1 u 1 1 ”b .j 


F 


(3.2.50) 


m -1 » i 

Bru(kT)z(Fi+m+.i-1 ) ] , 3=1, 2,..., F m^o 
^ B[u(kT)z(F¥m+3-l)3=o mf$J-U all 3 

=a^cr?,b5'”^ J m>1, all j 


F 
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Substituting (3.2.50) and (3.2.51) in (3.2.49) we obtain 

iiA(k,k+3i)= ^ [Vi®- 


11 N 


X T- S( z(N1^+i-1 )z(NK+m+ j-1 ) ) ] 

. i=1 0=1 




1 

’”^1 )2 .2 


K N 

5 , A 


B( z(N;^+i-»1 ) z{rrk+in+0“1 


(3.2.52) 


Bqn. (3.2.52) therefore is the recursive equation for 
evaluating R^(k ,k+m)| , Ih the limit, as N-J^oo 

RA(k»k+m)| o, it is seen from eqns. (3.2.46) and (3.2.52) 
that the ’autocorrelation function of the output of the coarse- 
sampled -data filter of the first order does not equal that of 
t?he output of the linear digital filter, even with infinite 
averaging time (N-*^ 

3.2.2 The Second Order Goars e~sarapled -data Filter; 

' # 

A linear second order digital filter is given by the 

difference equation; 

y(lT) = a^^u(lT)+b^y(TrTT)-t-b2 y(I^‘T) 

For convenience, let b^=o. Hence y(kT) is given by 

y(lT) = a^u(lT)+b^y(i::5'T) 


(3.2.53) 
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Kirlin’s coarse sampled data filter replacing e linear 
second order digital filter is given by the eqixations. 

Z(Nl+k)=sgn{a^u(ri)+n(m+k)+b22(Nl+k-2) ■ (5.2.54a) 

y(lT) = -^ ^ z(Hl+k) (3.2.54b) 

k=1 

The input u(lT) is as defined in eqn. (5.1.14). The period 
of the dither n(Nl+k) is T’ seconds, n(lTl+k) is uniform3.y 
diacributed in the range (-1,+1) so that 


a«u(lT)_^ + Ib^i 4 1 


O' ' 'max 

Sb 2(>-1 
Pn(n) =5 - 4 ' 


(3.2.55) 


(3.2.56) 


Fig. 3.2.2 shows a schematic diagram of the second order 
coarse~sampled~data filter. 


(a) The . hefcerrainistic Input : The deterministic input 

u(lT) is such that the eqn. (5.2,55) holds. The output of the 

coarse sampled data filter is given by eqh.s (3.2.54) 

(i) The probability of z(ITl4k) ; The initial values of 

z(llHk~l) and z(Nl+k-2) at l=o and k=1 are set such 
that the initial stored values z(o), z(-l) are each 
+ 1 vrith equal probability and are mutually 
independent i.e. 
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B{z(Nl+kJ=B[agn(a^u(lT)+n(Nl+k)+b 22 (Kl+k-. 2 ))] 
=aQu(lI)+b2S[ z(Nl+k-2) ] 


l-b, 


k-M k~1 






♦ Ct _ * 


o' l-b. 


(3.2,62) 


^ u(l-rT') (b^^^ for odd k: even B . v(®') 

r=t 

i£±l VI . - N+1 

1 -b^ 2 ^ l-b, 2 

=aQu(lT). + b2 ‘ .a^ * 


1 , 

Z, u(TIrT)(b„ 2 for odd k; odd N'. 

r=1 ^ 


(h) 


1 -b !*■’ , 

T -b — ■ + V •%• -jzi 


Sfi 


7 a(I-rT) for even kjeven N(c) 
r=] 

■ k. , • N+1 

l-b, T' , 1 ~b, 2 

=a_u(lT). —T— ^ +b^''^ .a . — r~r 

0 '' 1-bo 2 o 1-b;k 


1 -r'l 

r u(l-::rT)(b 2 ~ ) ' for even k; odd N(d) 

r=1 

( 3 . 2 . 63 ) 

Bor the two cases, (i) N is odd (ii) N is even B(y(lT)) is given 
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S(y(lT))= f _ _E(z(Nl+k))+ B(z(Nl+k))]N odd (a) 


k=1,3,5...N 


lc—2 f f • • ♦ # W"* 1 


= H E(z(Nl+k))+ TL B(z(Nl+k))] N even (^>) 




k=2,4,...,F 


(3.2.64) 


Theref ore , 


B(y(lT)) = # [ 




N ■- 1 -b. 

4 , 

N+1 


.N- 






F+1 

(l+b^-ab^ 2" )+ 


1-b- 2 N+1 1 N-l ^ , 

a 1 (l+b^-2b, 2 ) T u(l~rT)(b„“r‘ ]odd N (a) 


N. 


1 


. ,a-u(lT) a b^uClT) vr/o 2 

1 [41 .N . -g 2 1 . 4 (t+bj(l-bp)+ . 

^ ^ "^2 (l-.bj2 ^ 2^2 o ( 1^)2 


. 1 lAizll 

( 1+bp) T_ u(l-rT)(b52 ] even N 

^ r=1 ^ 


(b) 


(3.2.65) 


In the limit as IT-vcofor infinite averaging £ (y( IT) ) tends to 

a u(1T) 

B(y(lT)) = - 4 '- ;— — for N either odd or even (3.2.66) 


N--><o 


Comparing the RES of eqn. (3.2.66) with the o/p of the 
corresponding second order linear digital filter (3.1*50), 
it is seen that S(y(lT)) is not close to the desired 

value y(lT) . 
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(ii) The Time Autocorrelation, of y(lT): for the second 
order coarse sampled data filter is defined in a 
manner similar to that for the first order filter, 

K 

Z y(M)y(E+5T) {5.2.67) 

, xC““0 

li„ 5 K 1 N N 

"" E^ooST Z1 "T S z(Nlc+i)Z(zk+m+ 3 ) 

(3.2.68) 

E(I^^(m))= ^ ^ B(z(Nk+i)z(Nki¥4-3)) 

(3.2.69) 

E^^ dr 21-2“ Z1 ^ B[sgn(a u(kT)+n(Nk+i) 


+ h^z(Nk+i-2) ) .sgn(aQu(k+mT)+n(Nk+m+j)+b 2 z(Kk+m+j- 2 ) ) ] 


(3.2.70) 

B(E^Cw)) can be expressed in the form of recursive equations. 


For m=o 


B(BA(o))= ^4^ [N+(B‘=^-N)a^u'='(kT)+a^b 2 u(kT)(N~ 1 ) 

w . lr:im TT 


E 


1 


2 rr\^ 2, ,2 , 


k=o N 


I 


N N . , N N 

f E(Z(Nk+i-2))+ f E(Z(Nk+j-2))4 +b; Z t, 

i=1 3=1 * 1^1 j=1 

B(Z(Nk+3-2)Z(Hk+i-2))] 


(3.2.71) 
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Depending on the nature of N, the appropriate two equations 
of equation (3.2.63) may be substituted in eqn. (3.2.71), 
to obtain E(R(l(in)) by expanding the last term of -eqn. (3.2*71) 
recursively* 

For mj^o : From eqn. (3.2.70) for m^o, 

K , IT N 

7 ^ / L E[(a u(kT)+b.z(nI':+i-2)) 

k=o i^ 3 = 1, ° 

(a^u ( k+mT )+b2 z ( Nk+m+ 3 -2 ) ) ] 

K 

^ [N^a^u(kT)u(k+mT)+N a^b- u(fe+mT) 

k=o IT ° ° 

H IT 

■ ’ T E(z(ITk+i-2))+u(kT) J B( s(lTK+m+3-2-) )+ 
i^ 3=1 

- ^ IT IT 

bp T 7" E(z(E’1i+i-2)z(Nk+m+3-2) ) ] (3.2.72) 

• i=1 3=1 

As IT is odd or even, the appropriate substitutions may be 
made from (3.2*63) in ecpi . (3.2.72) and the last term in 
eqn* (3.2.72) expanded recursively to give B(R,^(m)). 

v 

(b) The Random Input : to the second order coarse-sampled- 
data filter u(lT) is white, stationary and independent as 
described in eqn. (3.2.14). In addition, its range is such 
that eqn. (3.2.55) is satisfied. The statistics of the 
output y(lT) of this filter are derived thus: 


lim 1 




lim 1 
K+1 
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(i) The conditional probability of z(Nl+E) ; The values 
of z(IJl+k-Ot z{Nl+k-2) are initialized at l=o, k=1 such 
tiiat z(o),z(-l) are +1 with equal probability i,e. eqn, 
(3.2.57) holds. Eqn. (3.2.54a) gives z(ill+k), a function of 
u(lT), n(]!'Tl+k) and z(]Tl+k-2). Hence the probability of 
z(Nl+k) conditioned on both u(lT) and z(Nl+k~2) is given by 


T) 

Z 


iz 

Nl+k ^Hl+k-.2,utlT) 


l+a u(lT)+b„z(Hl+k~2) 
(1 j z(Nl+k-2),u(lT) )= 2 ^ 


'Nl+k ^Nl+k 


(-lj 

-2,u(l 


(3.2.73a) 

1-a^u(lT)-b„z(Hl+k-.2) 

z(Nl+k-.2),u(lT))= — 2 2 

T) 2 

(3.2.73b) 


The unconditioned probability of z(Nl'-k) may be obtained from 
its expected value, evaluated next. 

(ii) The Expectation of '}(1T)*. Using eqn. (3.2.54) 

E(J(1T) = X B(z(Nl+k)) (3.2.74) 

^ k=1 


B( z(Nl+k) )=B[Bgn(a^u(lT)+n(Nl+k)+b2 Z(Nl+k-2))] 

= ii^+b 2 B[z(Nl+k-2)] (3.2.75) 


The result is as given in eqn. (3.2.63), u(lT) and u(i-rT) 
being replaced by the substitution made as before:- 
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Tj a N a ts 


IT4- 1 




E(y(lT))= ^ 


n“? 5 1 ”*^.5 2 

P ! A 2 \ _ 2 


2 (l-b^) 


F+1 


•2-(l+bp-2bp • )+a . ^ 

^ ° (l-bj)^ 


1 lK^(r-1) 

(1+b -2b ) I b_ ^ ] odd 

r=l ^ 


aQb2 


^2 (l-b2)' 




1-b, 


2 


w/o 1 N '(r-1 ) 

(l4-b2)(l-b//2) 2: b..— 

r=1 


(l-b2)‘^ 

] even N 

(5.2.76) 


E(y(lT) = -2-Li 


0.2,11) 


If .u(lT) is zero mean, SHS of (3.2.77) is zero, tallying 
with E(y(kT)) obtained from equation (3.1.50). But if 
7^0 (3.2.77) is seen to be greater. 


Evaluation of probability of z(Nl+k) ? 
bet p (1 ) = P 

na+K 


E(z)= p+ (~l)(l->p) 

p ( 1 ) = l lMt J j . (3.2.78) 

‘“ei+k 

Since E(z) is already available from eqn. (3.2.75) and 
•• (3.2.63) » 'P' may be evaluated. 



(iii) The Ensemble Autocorrelation of y(lT) is given by 


Ha ( k , k+m } (y { kT )y ( k+ml ) ) 

1 1 M- 

“ "T / T S ( z ( -nc4 i ) z ( ITk+m+ j ) ( . 2 . / 9 ) 

n i=1 3 = 1* 

m-o 

, IT F 

RA(k,k) V J E( 2 (Fk+i)z(Kk+ 3 )) (3.2.80) 


. F N 

= [F+ T E(z(Nk+i)z(Fk+ 3 ) )] (3.2.81) 

ir i=1 4^ 

iT^O 

1 N N 

= — *■ [f+ 3 7" Efsgn(a u{kT)+n(Fk-t-i)+b^z 

ir i=1 3=1 ° 

17^3 


(Fk+i~2) ) .Sgn(a^u(kT)+n(Nk+3)+b2a(Nk+3-2) ) ] 


^ [iI+(N''-F)a^a2+2a^b2(N-1 ) "L S(u(kT) 


F 


5 I'l IT 

z(Fk+i-2))+b^ 7 7 B(z(Nk+i-. 2 )z(nk+ 3 -.^) )] 




For odd N or even F, eqn. (3.2.81) reduces to 

a 

H^(k,k)= ^ [F+(I'I^-F)a2o2+2a^b2(lT-l) j^bf 


(3.2.62) 

2 

a 0 
c 2 


H-t- 1 
2 




(b,+b^--2b^ ^ ) +b^ H H B(z(KE+i-2)z(iTk4-3-2))] 

^ ^ i=1 3 = 1 

17^3 
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a 

o u 


•l5.[N+(N2-N)a2a2+2a b (N-1) (N-2)- 

0 u o 1 ^ 1-b2 (i-b^)2 


-bp)' 


.(l+bp) 




J I 

7 > E(z(Kk+i-2)z(Nk+3-2))3 for even K 

>1 

(3.2.83b) 


Thus it can be seen that with N'-^^sofor infinite averaging, 
' % 

eqn. (3.2.83) reduces to the recursive equation; 


lim 




U 


2a^b 

o 1 u 

1-b^ 


. lim -2 

N-^o<j 2 


N N 

z: z: 

L=1 3=1 
i?^3 


. . E(z(lTk4-i-.2) z(Nk+3-.2)) (3.2.84) 

For m^o eqn. (3-2.79) can be rewritten as, 

' i ^ I 

RA(ic,k+m)= T ^ B[(sgn(a u(kT)+n(Nk4-i)+b« z(Nk+i-2)). 

» ^ N i=1 3=1 ° 

(sgn(aQu(E+HT)+n(Nk+m+3)+b2 z(Nk+3-.2))] (3.2.85) 





S[ (a^u(kT)+b2z{Nk+i-2) (a^u(k+mT)+b2z(Nk+rn+;;- 2 ) } 


(3.2.86) 


= ^ [a b-N T B(u(k+mT)z(Nk+i-.2) )+a_b_N 
if ® ^ i-^1 o 


N 

H E(u(kT) 
0=1 


5 J N 

z(j!fk+m+ j-2)+b^ Z-. E(z(Hk+i->2)z(I'I&m'4-3-2) ) 3 

i=1 3=1 


(3.2.87) 



Thiu recursive equation may be evaluatei for 
and for different values of m. 


or even 


o 


1 


Suti.fTQ ry : The coarse-saii\pled-~data filter usin^fi: an RO avesrager, 
hes I een introduced by Kirlin [ 7 ] to replace a linear digital 
filter implementing by a particular difference equation. The 
expected value of its output y(lT) is said to be equal to 
that of the output of the linear digital filter. A saving 
in hardware requirements is made with the use of this 
non-linear filter. In this section (3. .2) the first and 
second order course sampled data filter, using digital 
counter averagers, are discussed, and the statistics of their 
outputs are compared with those of the outputs of the 
corresponding linear digital filters. It is seen that with 
this form. of the coarse sampled data filter, thf^* statistics 
of y(lT) do not approximate those of y(lT). 

3,3 The N-bit Feedback Coarse-Sampled-Data Filter: 

j.n the previous section, the coarse-sampled~data 
filter was discussed and it was shown tiiat it cannot be used 
to replace the linear digital filter, as such. A modification 
of Kirlin' s filter leads to the N-bit feedback coarse- 
sampled-data filter, the mean and covariance of the output 
of which can be made to match, as closely as desired, the 
mean and covariance of the output of the linear digital 



The tiller conf igaira tion is developel as 

follows : 

The oatput y(lT) of the coarse-sample ]- iata filter 
is obtains ’’ by lirst resetting the digital c'''‘'intL!’ aver'tgeT- 
to zero. Tile counter then upcounts or downcounis as 
z(lTl-t-k) is (tve) or (-ve) respectively, for 
Subsequently the sum stored in the counter is normrslized 
'by dividing by N, If it can be ensured that over 
the outputs of the comparator, z(Wl+k)'s, are Identically 
distributed ami tha.t E[z(Nl4-k) ]=E(y (IT) ) , then Mie expected 
value of y(lT) equals the expected *a.lue of y(lT). Th% s 
objective being satisfied, the respective au ^ ucorrelation 
fixnctiPn's of th*^ outputs of the two filters, th.--' linear and 
the non-linear, may be compared. 

To ensure identical distributions of z(nitk), k“1,...,jT, 
(the outputs of the comparator), the feedback for n filter of 
order M, contains, not fi elements as in the case of kit lln's 
coarse-sampled-data filter, but MW elements ' where W relates 
the dither rate T' and the input ^output sampling rate T’ 
(T*=T/W). The MW elements of the feedback, z(o , 
z(-MW4l) are set at l=o to +1 or -1 with equal probability 
and are mutually independent, tience for a linear digital 
filter of order M given by 



3 . 4*5 


y(lT)= 2 - a.u(l~iT) + ^ b.yd-jT) 


( 3 * V • t ) 


1=0 3=1 - 
the X-bit fecibac}' coarse-samplcd-data filter is given by 


the rel-'^tions 


M 


z(Wl+H.)=sgn( 71 a.u(l-iT)+n(lTl+k)+ 71. zCNl^j+k)) 

f=1 

(3. 3.. "’a) 


1=0 


A 1 

y{iT) = 1 7“ z(Ni+k) 


(3.3.2b) 


As in sec. 3.2, the dither, lying in the range (-1,+1), is 
such that 


n M 

[ 21 a.u(l-iT) ] + ■7_ 1 b.. j ^+“1 


and 


1=0 , 


M 


max i=i 


(3.3.3a) 


[ Z. a.udni) ] 

/ i=o min 


M 

S'**' 


^-1 


(3.3.3b) 


Choosing such a range for 'n* ensures that the input to the 
comparator, over a period T, is (+ve) as wv.ll as (-ve), 
avoiding saturation of the output of the cvxic r- *, which 
makes an averaging of the output of the comparator meaningless. 

The dither noise n(Fl+k) is independent between 
samples, is uncorrelated and is uniforaly distributed in 
the range (-1,1), i.e. if Pj^Cn) is the probability 
density of the dither, 



(3.3.4) 


p^(n^ = i -14 n <1 

Ihe output stativstics are e^'/^luatod for t!,.o f irat 
and second oi-dor filters, for the n^'ses of d ctera.ini at i c 
a nd ra nd o m inputs. 

3 . 3.1 xhe First Order .i-hit feedback Filter : 

The first order linear digitcsl filter follows the 
difference equation 

y(lT) ^ aQu(lT)fb^y(l-fT) (3.3.5) 

The corresponding N-hit feedback f liter (f i.-,. (3.3.1) 
is given by the equations 

z (Nl+k)=sgn(a^u(lT)+n(Wl+k)+b^ z(Ni~+k) (3.3 .6a ) 

A 1 

' y(w) = i Z (3.5.61.) 

k=1 

The added dither n(Nl+k) follows equations (3.3.4), 
with M=1 . 

(a) The Determinist- • Input: Ltt u(lT) be deterministic, 
its range obeying eqn. (3.3.3). 

(i) The Probability of z(lIH-k): The probability of z(Nl+k), 
conditioned on z(Nr-T+k) is 



he liuv^nr digital 
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3 .4<: 


1) 


i 1+a u(lT)+b 

(1 t z(: :-1+^))= 2 I 


-f ' “lil- 1-f k 


(•5.3.7) 


P 


L-i-k ! Vl-1+k 


i-a u(i:)--b.z(:':i~i-k) 

(-1! z{!Tiri+k))= — 2 — 1 


2 


(5.3.8) 


The unconditional probability p(z(Nl+k)) is given by the 

\ 

recursive equation 


z 


Nl+k 


( 1 ) = 


1+a u(lT)-b 


1 


^1 Pz 


( 1 ) 


Nl~1+k 

2^ ^ 


1=0 


1 


(3.3.9rO 


Similarly, 

P,' 


ITl+k 


1~a u(lT)-b. 

(„0= 2..^ L+ 


(- 1 ) 


Nl-1+k 


1 1-a u{iT)~b , . b 

_ ^ o 1 , ^ b"^”’^ + -4 


1+1 


1=0 


1 


(3.3.9b) 

Thus z(Nl+k), k=1,2,...N, are identically distributed, 
and independent of 'k' . 

(ii) The Expected Value of y(lT): 

Eqn. (3.3.6) gives the input-output relation for the 
first 01 x 3 or K~bit feedback coarse*sampled-data filter. 

Using eqn. (3.5.6b) 



(5.5.10) 


S(y(lT)) = 3 21 E[z(tll+k)] 

k=1 

From eqn. {3.3.6ba) 

B(z(Kl+k))=B[Sgn(a^u(lT)+n(lJl+k)+b^z(NirT+k)) ] (3.3.1 ! ) 

=E [a^u ( IT ) +b ^ z ( nH'T + k ) ] 

=a^u(lT)+b^E(z(BTIT+k)) (3.3.12) 

E(z.cm+k))=a ^ u(iT)b^-i ■ (5.5.15) 

i=o ' 

¥e see that E(z(lTl+k)) is independent of k. Substituting 
E(z(Wl+k))in (3.3.10) we obtain 

E(y(lT))=F ZL [a. ^ u{lT)t}-^ ] 

^ k=1 ° ^=0 ^ 

' 1 T • 

= a T u(iT)b;“^ (3.3.14) 

i=o 

Comparing eqn. (3.3.24) with eqn. (3.1.9) we see that the 
expected value of y(lT) equals y(lT), the output of the linear 
digital filter, for a deterministic input u(lT). Comparing 
eqn. (3.3.13) with (3.1.8), we see, too, that 


B(z(Nl+k)) = y(lT) 


(3.3.15) 



(iii) The Time Autocorrelation of y(lT); 


Lef the time autocorrelc tion function of y(lT) be 


denoted iy 


RA ( m ) 


1 in 1 






y(kT)y(I:+ffiT) 


’^A(m) is a function of the added dither n(Nl4-k) 


(3. ''.16) 

Hence 


B(RJ.(m))= U-j- X- B(y(kT)y(lt+mT)) 


K 


To evaluate E(y(kT) y(k4^T)) 


For iD=o ;• 


(3.3.17) 


B(y(kT)y(kT))= ^ 


J. X 

? 2- X-- S ( z(Mk+i 1 z(Nk+ j ) ) 

ir o=i 


1 J- 0 *' JEf 

^ [ X. B(z'^(Nk+i))+ XL^Hzd'T^+DzCr'X-^j))] 


N' 


2 


i=1 


i=1 3=1 


(3.3.18) 


B ( z ( Nk+i ) a( i'Ik+ 3 ) ) 



=B[sKn(a^u(kT)+n(2ik+-i)+b^ z( Ilk-k i) ) . 


Sga(a^u(kT)+n(I1k+ 3 )+b^ z(IIk-H -3 ) )‘] 

= u^(kT)+a^u(kT) .b^ [r;(z(nK-! ^ i)+ 

a(Nk3T+3) )]+t'^3(^‘(’'^“1+j)2( T+i) ) 



3.5t 


fi(y(kT)y{kT)) = Jy [N+(N^-.N)y^(kT)] (5.3.l9b) 


In the limit as 

l^ooMyikT))^ = y^(kT) (3.3.20) 


For mj^o. 

B(y(kT)?(EfST)= -U £ £. E( 2 (Nk+i)z(HETS+J>) (3.5.21a) 

i=1 3=1 


= c r. E(z(Nk+i)z(Nf+i+i))+ z. 

r i=i i=i 3=1 

ii^3 

E(z(Nk+i)z(l!rk+m+3) ] 

B(y(kT)y(EfmT) )| = [N(b“+a^y (kT) Z u(kTmH-i) 

'm>o ^ ° ' 

I 

4-(N^-N)y(kT)y(k+mT)] (3.3.21b) 

, . )_4 |m|-1 . 

E(y(kT)y(k+mT)) = -1^ [H(b + a y(S+5T) Z. b5-u(lo::i:T) 

* 1 B <0 ^ i=o * 

+ (M^- H)’ y(kT)y(Fi¥T)] (3.3.21c) 

Therefore, in the limit as N->oo, for infinite averaging and 

for all m E ( y ( kT ) y ( k+mT ) ) tends to 

E(y(kT)y(k+iT)) = y(kT)y(S+mT) (3-3.22) 



In tno 'lb I;— the expected vt? lue of P^(e') 

y 

( eqn . 3.3.17) tend 3 to 


lira 




lim 1 
K-^-^ K+1 


K 


k=o 


y(kl)y (k+mT) 


But by definition (eqn. 3.1.11), the RK3 of 
time autocorrelation function of y(lT). 


(3. ^-.23) is 


■? 1 


the 


i^co S(RA(m ) = Ry(m) (3.3.24) 

/ « 

This implies that if the filter be allowed to run for an 
infinite length of time to obtain each sample of the outt^ut 
y(lT), the expected value of the time autocorrelation 
function of its output y(lT) shall equal the time auto- 
correlation function of the output y(lT) of the linear 
digital fid ter . 

(b) The Ilandom Input ; Let u(lT), the random input to the 
filter, be defined as in eqns. (3-1.15). The input-output 
relationship of the i'irst order N-bit feedback coarse sampled 
data filter is given by the eqns. (3.3.2) and the range of 
u(lT) satisfies eqn. (3.3.3). The statistics of the output 
y(lT) may be evaluated thus: 

(i) The Probability of z(Nl+k): 2 (ili+-k) is given by eqn. 

a 

(3.3.6a). It is/.function of n(ni+k), u(lT) and z(Ll-l4-k) 





:fh' inrr<t a(li) is indej.>'~rient between ea'rples-, 

ru'nee u( I i ; r^nJ 2,( ’'l+r+G ) r^- c. , all s are ini epend ent - 

ihe pro ability ot z(ijl+’:), therefore, con'iticncl on 
u(l ) n ni z ( 1 -I > ) , is given by 

"i n ' ' b, •' ( ’ - 1 ’ ' 


P, 


'NH-hl "Wl-i+k ,^1T 


(1 jz(IIl-1+lO ,u(lT))^ 




z 


Nl+k 


m-l+k,*"!! 


(-1 


, 1 -a u{ lT)-b, z( ',1 -Hb ) 

z(Nl-.1+k),u(lT)= 2 ;_! 


The probability of z(Nl+k) conditioned on z(Nl-1+k) alone, ii 
therefore given by 


p(z(Nl+k) z(Nl-1+k))= 


p(z(Nl+k) z(Nl-H-k),u(lT)) . 


u 


1 


p(u(lT)). du(lT) 


( 3 , 3 . 26 ) 


( 1 z(Nl-1+k) ) 

^Nl+krwT~l+k I 


"Nl+k" 


(-1 z(NT-1+k)): 


1+aoP^+bi z(Fl~1+k) 


l-aohu"^! 




FI- 1 +k 


(3.3-r^Ya) 


('5.5-Y7b) 


where - E(u) ('5.3.27c) 

In the form of recursive equations, p(z(IIl+k)) is giv-ai 
•by 


Ik 



r/ 


Ua n -b 

) (1) = 2 ^ Pz 

^Nl+k ' %1-1+k 


^ , 1-i ^ ^ 

2 °I '*' 2 

1=0 


1+1 


Similarly, 


■“Nl+k 


(- 1 ) 




+ ’ . p 
1 z 


(- 1 ) 


Nl~1+k 




1+1 


^ ^ "g ‘ " 


1=0 


(3.3.2aa) 


(5.3.28b) 


(3.3.29a) 


(3.3.2913) 


Bqns. (3.3.28b) and (3.3.291 d) show that p(z(Hl+k)) is 
independent of k, and that z(Nl+k), k=l,...,N are identically 
distributed , 

Ai 

(ii) The expected value of y(lT): 

y(lT) is given by eqns. (3.3.6), u(lT) being rand or. « The 
expected value of y(lT) may be evaluated as follows. 


N 


E(y(lT)) = ^ 2lE[z(Nl+k)] 


(3.3.30) 


k=1 


From eqn. (3.3.6a) 

B(z(Nl+k))=B[Sgn(a^u(lT)+n(Nl+k)+b^z(IlTrr+k)] (3.3.31 ) 

= B[aQu(lT) + h^ z(NlM‘+k)3 
1 


ah >6, 
O 'U 1 


1-i 


(3.3.32) 



1 i pt; ' 


/ < ' ]t, , 


'"'/(r 1 : 




u j 

) 


'*„»i 2 i. 


l~i 


Comparing aqrs. {;<,3.34) ani ('’,1.16), wp* see f^hai. ;;)*■ 
of (3«'5.34) is equal to E{y(lT)). Ht^nou the expocto.i 
of y(lT) nnil y(lT) are equal, even if /-o. We I so 
o'''r:t,'f'vu that the- *U-Ib of eqn, (3.3.3?) is oqual to R(y(ri’;). 
The re fore , 

E(a(Ml+k)) M(y(lT)) (3.3,3a) 

(ill) The filnsemble Autocorrelation of ^(IT): io defined a a 


R^(k.,k+m) iify(k'J})y(k+mT) ] 


(3.3.36) 


h E i- ii(z(Hk+i)z(iiFH;+n) 

i=1 3=1 


?or Tii=o ;~ 

, N N 

RA(k,k)= V 21 T E(z(Nk+i)2(Nk+3) } 

y ft j:ri 


(3.3.37) 


(3.3.38) 


. K « N N 

•4- [ -r E(z‘^(Wk+i))t- 2- 21 B(z(fIk4-i)z(Nk+j))] 

if i=1 i=1 j=1 






(3.3.39) 


CENTRAL LIBRARY 
Acc, No. Ak Ew.TS.«L. 



It ia known thrit z(Nk+i), z{Hk+-j) i/j for all are 
identically distributed, hence the second term in ( 3 . ’^• 39 ) 
may be evaluated thus: 

N N „ 

j- Z. E[z(Nk+i)z(Nk+j)J=(r-.N)E[2(Nk+i)z(Nk+j)]^j^^ 
i/j , , 

(3.3*40) 


= (N^-N )E[sgn(a^u(kT)+n(Hk+i)+b^z.(Nk-.1+i)). 

I 

Sgn(a^u(kT)+n(Nk +3 )+b^ z(nF-1 + j ) ) (3.3.41 ) 

= (N^-.N)E[(a^u(kT)+b^z(HSi:!'fi))(aj^u(kT)+b^z(NCT+j))]jL^j 


(3.3.42) 

=(N^-.N)[aQcr^+aQb^. E|^u(kT)z(ITkrT+i)+u(kT)z(Nkrr+j)| 

+b^ B(z(NErr+j)z(NErr+i))]^^^ 

= (E^-N)[a^o^+b^ B(z(NFrT+i)z(NErr4-j))]^^^ 


(r~H) ajo^ 


(. -2k+2 ^ 
2-2 ^ 


1-b 


1 


2k+2 


• . (k,k) = ^ [N+(H2-H).a2o2. J 

In the limit, as N->«>for infinite averaging, 

■ 2k+2 


liiiii 


RA(k,k)= 


(3.3.43) 

(3.'3.44) 

(3.3.45) 


(3^3.46) 
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For m^o 

, V 1 r J- ’ N N 

R^(k,k+m)= [ 4lB(z(Nk+i)z(lTE¥+i))+ X 7L 

^ i/ i=1 i=1 j=i 


B(z(Nk+i) z(Nk+m+j ) ) ] 


(3.5.47) 


= [lTb“ +(N^~K)'b“.a^cr^ , i* — ] m *;> o 

ir ^ 1 0 u ^^^2 

1 (3.3.48a) 

1 (k+m)+2 

= ^ [B +(N^-N)b)“‘^ a^a^ 2 — « ] m< o 

N*^ ^ 1 o u ^^^2 

(3.3.48b) 


In the limit, as N-><k? for infinite averaging, (k,k4-m) 

f * 

tends to a limiting value: 


l^^oo R^(k,k4-m) = aVd-b^)-^ m = o 


(3.3.49a) 


= a^o^ bj“Wl-b^^^'^“^‘^^)(t-b^5"^ m< o 

(3.3.49b) 

_ , m /, v2k+2v /, * s. 

^■o^u ).(1-Dd ni> P 

(3.3.49c) 

Comparing eqns. (3.3.49) with eqns. (3.1.25)» we see that 


N^<x> P'^dc^ls^+m) = Ry(k,k+m) 


(3.3.50) 


I 





where R {k,k-t-m) is the ensemble autocorrelation function 

Of 

of the output of the linear digital filter. 

3*3.2 The Second Order IT-bit Feedback Filter ; 

A second order linear digital filter follows the 
difference equation 

y(kT)=aQu(kT)+b^y(E:TT) +\)^ y(k^T) (3*3.51) 

f % 

The corresponding N~bit feedback coarse sampled data filter 
is shown in Rig. (3*3.2). The relation between the input 
u(lT) and the output y(lT) is given by the equations: 

7,(]!iri+k)=sgn(aQu(lT)+n(Nl+k)+b^ z(Nl^r+k)+b2z(N'l-2+k) ) 

. . (3.3.52a) 

y(ir) .= w 21 z(Ni+k) (7.3.52b) 

^ k=1 

Since this is the coarse-sampled-data filter 
implementation of a linear filter of order two, therefore 
as stated earlier, the feedback consists of 2R elements. 

At 1=0, the feedback contains z(o), z(-1 ) , . .'.z(-2N+l ) each 
of which is +1 or -1 with equal probability i.e. 

P (1)=P (-l)=i r=o,.1,-2,..,,-2N+1 (3.3.53) 

The added dither n(Nl+k) follows eqns. (3.3.3) and.X3.3.4), 
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with M=2. For almplicity, let b^=o. Hence, the non-lineer 
filter is given by the equations 


Z(Nl+k)= 8 gn.(aQu(lT)+n(Nl+k)/b 2 z{NT^+k:)) . (3.3.54a) 


y(lT) = i 21 z(Nl+k) (3-5. 54b) 

k=1 

(a) The deterministic input : The input u(lT) is detcrmiuib s. i . 
its range obeying eqn, (3.3.3) for M-2. 

« 

(i) The probability of z(ia+k} ; Th^ probability of z(NHk), 
conditioned on z(Hl-2+k) is found to be 


. _ (lU(NHT+k))= — ^ 

^Nl+kj ^Nl-2+k \ 


Ha .u ( IT ) +b« z ( Njt'~^+k) 


*(3.3.55) 

- 1-.■^u(lT)-b^z(MT^+k) 

_ (-1 U(Hl-2+k))= 2 “ 

^Nl+kj 'Nl^+k \ 

' (3.3.56) 

The unr- nditional probability of v.«'" .Mf) is ^ « 

recursive equation 

Ha„u( iT)-b 


P. (1) 

^Nl+k 


+ b^ P (1) (3.3.57a) 

2 ^%ir?+k 


Ijz HaQu(I^iT)-b2.b2 ^ 1/2+ H 

^ / 

i=o 


for even 1 


(3.3.57b) 
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1-t 


X . l±i 

^ 2 5 (»>2 2 ) 


for odd 1 


Similarly, 
Pz (-1) 


(5.3.57c) 


'Nl+k 




1/2 


“ ? (’•^o“(5^iT)-b2)b| + i(b’^V2, 


1~1 

2 


2V-2 

for even 1 (3. 3. 58b) 

l±i 


it- ? ('-a<,“(l-2iT)-b2).bi + i; b^^ 


for odd 1 (3.3.58o) 

Thus z(Nl+k).k==i ? w 

U2»...,N are Identically distributed 


and 


indepwndoiit of 'k', 

/* 

(ii) Exp ected Value of rj^. ^ , 

— using eqn. (3.3.54b ) 


we obtain 

,A 


N 


E(y(lT)) = i ^ B(z(Nl+k)) 


k=1 


(3.3.59) 


Using eqns. (3.3.57) and (3.3.58), 


1/2 


E(z(Nl+k))=aQ u(r:?iT).bj for even 1 


(3 f 3.60a) 



3.62 


(l-t>/2 . 

=a X. u(l--2iT).bp for odd 1 (3.3.60b) 

i=o 


Therefore, 


E(y(lT))= ^(N).a^ 


1/2 . 

;2l ’u(!C::nT) .ht 

i=o 


for even 1 

(3.3.6la) 


( 1 - 1 )/2 ^ ^ 

= Z- b| 


odd 1 (3.3.6lb) 


= E(z(Nl+k)^. ‘ (3.3.6lc) 

Comparing eqn. (3.3.61) with the RHS of eqn. (3.1.50) we 
see that the expected value of y(lT) eqmls the output of 
the linear filter y(lT), when u(lT) is deterministic i.e. 

. ■E(y(lT)) = y(lT) (3.5.62) 

Since E(y(lT)) = B(z(Nl+k)), therefore when u(lT) is 
deterministic 

B(z(Nl+k))= y(lT) (3.3.63) 

A 

(iii) The Time Autocorrelation Function of y(lT); is defined 
as follows: 

Hj(m) = ^ £ ?(kT)y(ET5T) (3.3.64) 

BA(m) is a function of the added dither n(Nl+k). Thus 
B(R5(m))= ^ B(y(kT)y(S5¥T) ) 


(3.3.65) 
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The Inner term is given by 

N i=t pi 


)z(NFi^+4)] 


^ N E 

? 5, ^ ®(!!(Nk+i)z(lffi5-+j)) 

^ J I 


For m=o 


(3.5*66) 


(3.3.67) 


^ 1=1 d=i 


z(Nk+o)) 


(3.3.68) 


it? ^ (Nk+i))+ g- ^^E(z{Nk+i)z(Nk+j))J 

(3.3.69) 

" lf2 ^ ^,®lf«“K“('=f)+“('nc+i)+b2z(NE:?+l)). 
i?^d 

Sgn(a^u(kT)+n(Hk+j)+b2z{BK5+j)) ] 

+b2z(NF:2-fd)) j 

' ? jSi 

i/d . 

+*(N^+3))»-b^ E(z(lfkr2+i)z(iiSI5‘+j)) j 



J.64 


1 r I- J- 9 5 

- J- [Nt r 21 


JJ? *0^^ (kT)+2aQu(kT) .b^yCEI^T)^ b^* 




E(z(NE::2+i)2(NE::?+:j)) ] (3.3.70) 


= 1 
N 


J [N +(N^-N) y^(kT)] 


Therefore, 


K 




(3.3.7t) 


fN+(N^-F)y^(kT)] (3.3.72) 


In the limit, as N-^o®, B(R^(m)) tends to 

V 

N^ooE(R^(o))= i y2(kT) 

k'^'O 


(3.3.73) 


But the RHS of eqn. (3.3.73) is the time autocorrelation 
R (m) of y(kT) at m=o 

V , • 


N^oo = »y(°) 


(3.3.74) 


For m^oi- 

Zy/H m V A 


N N 


B(y(kT)y(k-fmT) )=B[-45- ^ 21 Z(Srk-f i) z(Nk4-m'f j ) ] 

i=l 3=1 


[ X B(z(Nk+i)z(l!lS+m+i))+(N^-N)B(z(Nk+i) 

N i=1 

Z(Nk+m+3)) w3 (3.3.75) 

't m/2 - 

• ; . ' +y(kT). XI ^2 
* q=o 

+(l!l^.-S)y(HP)*y(E7nrT)3 m>o and even 

(3.3.761*' 


N 


1 
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^ [H(bl“l/2 +y(ETmT) b| u(Er?,T))^.(H2_N). 


q=o 


y(kT)y(S+ST)] in<o and even 


{3.3.76b) 


If ffl la odd, z{Nk+i), z(ITk+m+i) are always independent, for 
all i , j , therefore 


H N 


l(y(k!r)y(k+mT))= ^ ^ ^ E(z(Nk+i)) .E(z(hS^^-4)) 

IT i=l j=i 

(3.3.77a) 

“ ^ ‘ .y(kT) .y(k+m!r)=y(kT)y(k+inT) 

(3.3.77b) 

Using eqna. (3.3*76b) and (3.3.76), we see that, for 
all m, in the limit, as N— > otf 


11m 




(3-3.78) 


lim 


^ y(W)y(BST) all . 
k*o 


K 


= V“) 


<3.3.79) . 


where By(m) is the time autoooi- - function of y(k!E) . 

(b) The random input : The input u(iT) ie zai^om. Its range 
obeying eqn. (3.3.3) for M=2. 



(i) The Probability of 2 ;(Nl+k): The probability of z(Nl+k), 
conditioned on both u(lT) ani z(Nn?+k) is found to b# 


'Nl+k 


, i , 1+a u/lT)+b5z(!rr::5^<-k) 

2 (lU(Nl~2+k),u(lT))= 2 2 1 

^ia-.2+k,«lT 2 

(3.3.80) 

J „ 1^^u(lT)-b.z(NrT+k) 

^ (-1 z(HX:2+k),u(lT))= 2 2 

f m I i. 


^Nl+kl *^Nl-.24-k*“lT 


p{z(Nl+k) 


(3.3.81) 

z(Kr:5 k))= j PV... 'i -fnir^'+k), u(tT;)- 

u 

.p(u{lT)).du(lT) (3.3.82) 


Therefore, 


■P, 


"Nl+k 


"Kl-2+k 


/ 1 , _ 1+a n +b,z(Nl-2+k) 

(ljz(NlI7+k))= 2-2 .. 2 (3.3.83) 


, , 1-a n -b.z(NT:::?+k) 

Iz -rl- \2{Ni:::2+k))r- o u 2 
^Nl+k\^Kl^+k 


(3.3.84) 


The unconditional probability of z(Nl+k) is given by the 
recursive equation 

1+a«n„--b, 


Pg. ( 1 ) 

^1+k 


o'*u “2 . , ^ . V 

^ \ f J 

^ ^l^+k 




(3.5.85a) 


^2 ^ for even 1 

(3.3.85b) 


i=sO 


2 ^ ^2 1*^2 

i=*o 

for odd 1 


(3.3.85c) 



Similarly, 


z 


'Nl+k 


(- 1 ) 


1-a n^-b 

+ >=2 P,. ^ (-') 

Nl-2+k 


(3.?. 86a) 


'-®oV-’>2 Hi ,i .1/2+1 

2 ^ 21 ^2 even 1 

i=o 


(3.3.86b) 


1..a T) -b (^l)/2 i 

' — 21 b|+ 

x =0 ^ 


,( 1+1 )/2 


Odd 1 (3.3.86c) 

Hence z(Nl+k;) ,k=l ,2 , . . . ,]!T are identically distributed 
and ind ep eiid ent of ' k ' . 


(ii) The, Expected Value of y(lT) ; 

^ 1 N 

B(y(lT))= 22 E(z(Nl+k)) 


k=1 


(3.3.87) 


1/2 , 

E(z(lTl+k))= a^Ti^ ^ for even 1 (3.3.88a) 


x=o 

( 1 - 1 )/2 . 

k=o 


^^O^u 1?1'' "^2 ^ (3.3.886) 


Hence 


A 1/2 . 

E(y(lT))= a^T)^ '^^2 for even 1 

i=o 


(3-3.89a) 



3.68 


Comparing eqns. (3.3.89) e ni (3.3.88) with (3.1.50)»we 
see th^jt B(y(lT)) equals the expected value of y(lT), and 
that B(z(Nl+k))al 80 equals E(y(lT)'k i.e. - 


E(y(lT))=B(z(Nl+k))=E(y(lT) ) k=1,2,...,N (3.3.90) 

(iii) The Ensemble Autocorrelation of y(lT): is define! by 


jRA(k,k+m)=E[y(kT)y(k+mT) ] 

. N N 

=S[-V ^ H z(Nk+i)z(!TE+S+j)] (3.3.91) 

jr i=i 0=1 


11 ^ N 

xqA(k,k+m)j = -1^ [N+ Y. Z- ' Oz(Fk+j)] 

^ •ii=o A\f i =1 0=1 

' 

= [N+dT^-N) B(sgn(a u(kT)+n(Ek+i)+b„z(N^+i)). 

, . o ^ 

sgn (a^u(kT)+n(Fk+ j )+b^ 

z(i^+j)) 3 ^ 

±f^i 


= -^ [N+(N^-IT) a^E(u^(kT))+2a^E(u(kT)). 

b,B(y(E^T))+b^ B( z(irF^+i) z(nS^+0 ) ) ] 


= -U [N+(k^^-.E)B(y^(kT))3 (3.3.92) 

. N 

In the limit as tio 


= Ry(k,k) 


(3.3.93) 





Is Odd, and z(i^k+i) and y(kT) are given by eqns* 

(5. .5. 54a) and (3.1.50) respective!,, then z(irk+i) , z(if5T£r<+ j ) 
are mutually independent for all 1 , 3 . Similarly yCkT), 
y(k-iinT) are independent. Therefore using (3.5.91)# 

. N N ' 

fiA(k,k+m)= -ly 77 27 E(z(Nl:-i-l)).B(z(Nk+m+3) ) ( 3 . 3 . 94 ) 

E"" 1^1 3=1 


E(y(kT)).B(y(1?f5T)) 


o ( = o) for odd m 


(3.3.95) 


When m is oven , and k=2n or k=2n+1 and k4-m=2r or k+m= 2 r +1 

. il N N. 

RMk,k+m)= ~p [ X E(z(Nk+i)z(Nk+in+i) )+ X "21 
^ i=1 • ' i=1 3=1 


B( z(lTk+i) z(Nk+m+ 3 ) ) ] 


(3.3.96) 




N' 




z: z: 

i=1 3=1 
ipi 


E(z(lTk+i) .z(Hk+ 3 ))] 


(3.3.97) 


RA(k,k+m) = [N.b^®/^+ b2' ^ ^ . (iT^-N)B(y^ (kT) ) ] 

^ I m > o N " 

(3.3.98a) 


Substituting from (3.1.54a) in (3.5.98a), we obtain 

2 n +2 




{k,k+m)| = ^ 

* m ^ o F 


1-b, 


1 -b: 


1 


(3.3.98b) 
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-.,k4-m) 


m< o N “ 1=1 j=i 

i,^o 

B ( z ( Nk+Ei+ i) z(Nk+m+ j ) ) ] 


( 3 . 3 . 99 ) 


R 4 k , k+K ) 


3, 1 

ia< o 




(3.3.100a) 

Substituting from (3.1.54a) in (3. 5. '100a), we obtain, 


RA(k,k+m)|= 


[N.bl,” 


1/2 


+(N ~N) .b 


fmj/2 


- ^2r*f2 

1-b^ 


(3.3.100b) 


Using eqns. (3.3.95) , (3. 3-.98b) , (3.3.100b), in the 
limit, as N — > 


R^(k,k+m) = Ry(k,k+m) all m. 

BUI4MMY: 

In this chapter the coarse-sampled-data filter and 

the h-bit feedback coarse-sampled-data filter have been 

investigated. It is found that the performance of 1he 

is not 

coarse-sampled-data filter^close to the performance of 

the output of the corresponding linear digital filter whereas 
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for both the cases of the deterministic and the random 

m 

inputs, the expectei value of the output y of the K-bit feed 
back sampled data filter is equal to the expected value of 
the output y of the linear filter. It is also shown that 
as i'J increases, i.e. as the dither rate T' increases for a 
fixed input sampling rate T(T'=T/lO» the auto-correlation 
function of y^the output of the N-bit feedback sampled -data- 
f liter asymptotically approaches the auto-correlation 
function of y the output of the linear digital filter i.e. 

BCy) = ii(y) 

and RK(k,k-Ha) = Hy(k,k+m). 

At steady state, i.e. with k tending to infinity, y{ll) 

1 T m 

becomes stationary i.e. (k ,k-Hn)=Ry(m) , therefore 

when thq N-bit feedback filter is allowed to run for some 
time, its output must become stationary, i.e. when the 
non-linear filter reaches steady state, the power spectrum 
of its output may be defined. Because the autocorrelation 
functions of the outputs of the linear and nonlinear filters 
are equal in steady state,- therefore the power spectra 
must be equal also i.e. 

S* (CO) = Sy(ljO ) 

When a comparison is made of the hardware requirements for 
the coarse-sampled-data filter and the N-bit feedback 
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that at time tH41+k, k=1,..W , the values 5i(IJs-ih')» 

3-1-1, 1-2, 1-M are stored in tiic fi^edback register, i.e. 
the* feedback register length must be M. It can be shown 
that the expectation of the oucput of this filter equals 
th‘- t Oj t-he output of the linear digital filter. However, 
the autocorrelation of its outi^ut is larger then tb-t of 
the output of the d-bit feedback filter, while being lesser 
than that of the output of thi coars'’-sampled-data filtei'. 
Here, a tradeoff between the cost of the filter and the 
degree of the accuracy of the output* of the filter, must be 
made, in deciding on the nonlinear ff^ter. A st-ady by 
simulation has been made only of the h rlt feedback sampled- 
data filter, however. 


/ 



CHAPTER 4 


RESULTS OF SIMULATION OF N-BIT FEEDBACK 
COARSE-SAI^LED-DATA FILTERS 

In Chapter 3, a study was made of what was celled 
the N-bit feedback coarse-sampled-data filter which was 
developed to replace a linear digital filter. The advantage 
was a simplification in hardware. In theory (Chap. 3.3) 
it was found that the output of this nonlinear filter 

I 

could be made to approximate as closely as desired the 
output of the corresponding linear digital filter in the 
correlation-invf riance sense. 

This chapter outlines the study made of the 
performance of the nonlinear filter by simulating it on 
the computer DEC-lO, which is a 36--bit wordlength machine. 
Sec. 4.1 discusses the program written for the DEC-10 
to simulate the nonlinear filter. Sec. 4.2 discusses the 
performance of 1st and 2nd order nonlinear filters as 
obtained by simulating them on the machine. 

4.1 Simulation in Fortran-IO; 

A program NLP (Non-Linear Filter) was written in 
the computer language PORTRAN-IO to simulate the 
nonlinear N-bit feed-back coarse-sampled-data filter on 
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the DBC~10 computer, A listing of the program is 
provifiel in Appendix A, The variable names used for 

m 

some pertinent quantities, input, output and intermediate, 
are explained as follows: 

The linear difference equation from which the 
nonlinear filter is derived is given by 

y(lT)= u(l-i+fT)+ b. y(l- 3 T) (4.1,1) 

i=1 ^ 3=1 

Consequently, the nonlinear filter is given by the equations 

RU RY 

Zhat = a .u(l”i+i T)+ n(Nl+k)+ b .zCNl-j+k) 

i=1 3^ ^ 

(4.1.2) 

z(Nl+k) = Zhat/jZhatj (4.1.3) 

• That (IT) = ( jZ z(lll+k)) LTMN/N (4.1.4) 

k=1 

I" 

Here n(Nl+k) is given by 

n(Nl+k)=(lIMN) . (RFB(lTl+k)-0.5) .2.0 (4.1.5) 

where RFB(Nl+k) is a pseudo-random number lying in the . 
range (0,1) and generated by the computer library prog3?am 
1G-UB. The factor ’2.0' normalizes the value (RFE(Nl+k)-0.5) 
to the range (-1,1). LIMN is the magnitude of the 


I 
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maximum excursion of h(Nl+k) either above or below zero. 
In terms of the filter coefficients and the input u(lT) 


\ 


RU RY 

LIM > [ 21. u(l-i+fT)] + ^ 

i“1 rnctX J = 1 



(4.1 .6a) 


-LIMN ^ [ 


RU 

21 a. u(l-i+lT)] . 



(4.1 ,6b) 


Other variables used are T,N,MU and M. They are explained 

thus : 

T : is the input/output sampling period. 

N : relates the period T' of the dither n(Nl+k) to the period 
T of the input, thus T'=T/N. If the order of the linear 
digital filter is RY, the ntunber of elements in the 
feedback register is (RY)(N). 

MU: (MU),T is the time allowed to elapse after which the 

filter is assumed to be in steady state i.e. the output 
of the nonlinear filter is assumed to be stationary. 


M : At time (MU).T, when the filter is assumed to be well 
into steady state, the previous M values of say, the 
output i.e. Yhat(MU-M+1T) , ..., Yhat(MU) are used to 
obtain the power spectrum of the variable, which in 
this case is Yhat. 



Before the first value of u(lT) ie input to the eystear, 
the feedback (z( j) , j=o,-1 , -N(RY-l)+1) is initialized. A 
pseudo-random number YY, generated using the external 
defined function idljDYl(}-) on the computer, is used (YY=RirDY1 (X) , 
RIJDYI(X), and hence YY, lies in the range (0,1), To 
satisfy the condition 


P,. (1) =p,,(-o 

0 3 


3= 0,-1,...,-N(lY-1 )+l) 


(4.1.7) 


z( 3) is initialised to +l(-l) as YY is greater than or 
equal to (less than) 0.5, or vice versa. The filter may 
now be allowed to run, generating an output value Yhat(lT) 
for every input value u(lT). It is seen that, to reduce 
the possibility of any correlation between the dithering 
noise n(Nl+k) and the initial values of the feedback 

z., o=0,-1 , ...,~N(Ry-l)+1 , their generators in the system 

3 

must be two physically different, independent random 
number generaters, as GGUB and BNDY1 are. 

The M-point power spectrum of a particular variable 
is obtained by using the subroutine MBil, which uses the 
Maximum Entropy Method for the derivation of the power 
spectrum of a variable from cdrrupted data. 
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4.2 Simulation Results; 

The computer program used to test the performence 
of the nonlinear N-bit feedback aampled-data filter has 
been described in the previous section. The program was 
used to simulate first and second order nonlinear filters. 

The First Order Pilter : follows the linear difference 
equation 

y(kT) = u(kT) + ay(k~lT) (4.2.1) 

where, to ensure the stability of the filter, | a \ . 

The performance of the first ord,er nonlinear filter was 
tested for various inputs. 

(a) Input Zero: When the input to the nonlinear filter is 
zero, the output y(lT) is given by the eqtaations 

RY 

z(ill+k)= Sgn(n(Nl+k)+ 21b. z(irri.+k) ) (4.2.2) 

i=1 ^ 

' f\ , R 

y(lT) = 1 ;^z(Nl+k) (4.2.3) 

" k=l 

The output ’therefore shall be the natural or the 
initial condition response of the nonlinear filter. It is 
zero mean and converges to zero with time because circulation 
through the filter causes the feedback to be weighted by 
which is a decaying term. This is illustrated in 
Pig. 4.1 » where the output was obtained by setting u(nT ) 




Fi g . 4 . 1 Initial cond i i * 'r r ^ 
first orde? N' v'i’S - ■ 
sampi«d-<!a/a hi-,- 
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to zero and allowing the filter to run. The output y(nT) 
ia aeen to be zeroinean» and close to zero. 

(b) Step Input: When the input u(nT) to the linear filter 
is a step function, the output y(nT) must ultimately 
reach a constant value given by 

y(nT) = (4.2.4) 

Because B(y(nT)) = B(y(nT)) 

Therefore B(y(nT) = (4.2.5) 

% 

For a particular u(nT) and a, the range of the dither 
LIMN is determined from 


LIMN > u_-_ + 

• v' max 

I“1 

(4.2.6) 

-LIMN ^ u„,„ •- 
. ' ^ min 

hi 


this case u„„^ = 

u„.^ = u^ 


max 

min 0 


LIMN ^ u^ + 1 

,“l 

(4.2.7) 

-LIMNZ u„ - j 
^ 0 1 

!“! 


LIMN r:- max [ 

Iv l“l| - K - 1 

(4.2.8) 


The maximum value that y(nT) may attain is 

K— I 

= i . N LIMN = LIMN (4.2.9> 
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Therefore If LIMN y- g “• , the output of the 
nonline^-.r filter shall saturate at +Lir4N. This implies 
th'^t the range of the input u(nT) is constrained by the 
inequality 

u 

max f u^ + a , u^ - a ]> (4.2.10) 

This is illustrated in Fig, 4.2 for which a=0.75. The 
output corresponding to an input =0.5 to the linear digital 
filter i.e. 2.0 LIMN=1.25. The nonlinear filter output 

saturates at 1.25. 

(c) The sinusoidal Input: The performance of the first 
order N-bit feedback sampled-4ata filter is evaluated when 
the input is sinusoidal of the form 

u(kT) = cos ) (4.2.11) 

the period of the input being AT seconds. Eqn. (3.3.lo) 
for the expected value of the output of a linear filter is 

i:(y(lT)) = b^ u(rriT) (4.2.12) 

i=o* ' 

The output y(lT) is therefore the weighted sum of siriuuoidal 
terms, each of period AT secs. Hence y(lT) must be 
periodic with period AT secs. 



y (nT) 
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Fig. 4.2 


Firtt ordar filter step response 
LIMN « t.25 , N . St2 


4,10 


The output of the N-bit feedback sampled-data 
filter, corrviBpond ing to the linear difference equation. 

y(kT) = u(kT)+ 0.75 .y(CTT) ( 4 . 2 . 13 ) 

where the porioi of the input is 16T, is 'iopictei in 
Fi:. 4.3(a). The output y(lT) is indeed seen to be 
sinusoidal with period l6T.Fig. 4.3(b) shows the power 
spectra of th s of the linear and non-linear filters 

and the povrer spectrum of the sinusoidal input, which has 
a period of 16T secs. The figure (4.3b) shows a peak for 
the same frequency l/l6T for the outputs o'f the linear and 
the IT-bit feedback sampled-data filters and for the input 
which is sinuso4.a..i.l . There is a finite dispersion about the 
peak frequency u'. each case. This is due to the fact that 
a finite number of points (M=80) were used to determine 
the power spdotram -in each ca.se. 

(d) V/hite Noise Input; Sec. (3.3.1) discussed in theory the 
performance of a first order nonlinear digital filter, to a 
TiJ-hlte noise input. Here we discuss the performance of first 
order nonlinear iiiiers corresponding to linear digital 
..filters ' given uj o^qn. v4.2.l) when the input is white noise , 

When simulating a first order nonlinear filter, for 
particular run of the simulator, the parameters kid and N, 



y(nT) 



Power spectrum 


4.12 



Frequency (rad/sec) 


order filters , output 
power spectra for sinusoidal input • 
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which have been previously defined in Sec. 4,t, have 
particular values. 

Fig. (4.4a), (4.4b), (4.4c) compare the power 
spectra of the output s of nonlinear N-bit feedback - 
coarse-sampled-ilata filters with those of the corresponding 
linear digital filters. For a=0.25, (0.5), (0.75) the 
values of 14U and W that give the closest pov/er spectrum 
matching are MU=512, (1024), (512) and N=5 12, (51 2) , ( 1024) . 

The parameter N relates the dither period T' to the 
input sampled period T (T‘=T/N). Increasing (decreasing) 

N implies that the dither period T' decreases (increases) 
for a fixed input sampling rate T. It is seen that by 
obtaining the output power spectra for a nonlinear filter, 
for a range of values of N there is a particular values of 
N that yields the most close matching of the o/p power 
spectrum with that of ;.he o/p of the linear digital filter. 

An intuitive explanation may be given. In +heory an 
infinite time i.e. N= co is required for perfect matching 
of the power spectra of the outputs of the nonlinear filter 
arid the linear digital filter. However, this is not 
practicable. Hence for low values of N, there shall be a 
finite error between the spectra. As N incre^Jaes, this 
error may decrease. An inherent feature in any finite 
wordlength machine is the introduction of round off /trim cat ion 



Power spectrum 


0.6 


0.4 


0.0 





Power spectrum 
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Fig. 4. 4(19) 


First order fitters » output 
power spectra . 



Power spectrum 
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y<kT)*u(kT) + 075 y (k-IT) 
MUs 512 , N =10 24 
LIMN »1.00 , T * 0.1 



Fig. A *4(0) First order fillers, output 

power spectra . 
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noiFie into the filter. As N Increases boundlessly, the 
eunulative effect of the cut-off/ truncation noise may 
inc^rease anl cause the error between the spectra of the 
linear an.] tlie nonlinear filters to increase. However, tiiis 
explanation remains an intuitive onr, as no sto-y, oitver 
theoretical or by simulation has been male to lack it. 

The initial condition response lies v’own as . 

Hence MU must be large enough for the initial condition 
response to be negligible. Therefore MU must be grea tea’ 

as Hoxirever ttie simulation results show optimal 

I ) 

values of MU are MU=512 for a = 0.25 and 0.75 and MU=1024 
for a=0.5 for each of which 'r. 0. 

Secontl Order Filters ; Pigs 4.5(a ) , (b) , (c) dp^dct the spectra 
for the linear digital filters corresponding to 

y''kT) = u(k)4a^y(i>TT) + a^yCt^T) (4.2.14) 

and the N-bit feedback sampled-data filters derived 1 rom 
them, both for white noise input. 

The simulation program is run to one to 

study the effect of MU and TI on tte performcncf' of the 

of the nonlin-iar f jlttu’s' ]M-r' orr.' nee 


filter. Tlie belr viour 



Power spectrum 
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Fig. 4.5(a) Second order fitters, output 
power spectra . 



Power spectrum 
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Fig. 4 .5(t)) Second order fitters output 
power spectra . 



Power spectrum 
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Fig .4 • 5(c) Second order filters , output 
power spectra . 



with N was consistent with the explanation given for the 
’'white noise input to first order filters” case. For 
a fixed MU, varying N over a particular ran^e, the 
difference hotv;een the povrer spectra of the linear and 
nonlinear filters is mininium for some II, above and lelo 
which the error increases- ?^g, 4.5(h) depicts the 
output power spectra for two different valut>s of 4, 
fixed MU. 


i‘o r 



GIIAPTEH 5 
COHO LtJ SIGNS 

In this chapter the results obtained in the preceding 
study are stated in brief and suggest iona for further work 
on the problem are mentioned. 

Lineer digital filters are characterized by constant 
coefficient difference equatiors. They, can be realized 
using the appropriate units of ■ digital hardware available. 
The coarse-sampled-data filter vrais introduced with the 
ifiea of re] lacing the linear digital.f liter with it. The 
expected value of the output of ‘the linear digital filter 
was said to be equal to the expected value of the output 
of the proposed filter. The advantage was that lees 
complex ha.rdware was required than in the linear digital 
filter case. However it was noted that due to the 
requirement of a high frequency dithering signal by the 
nonlinear filter, the input signals were constrained to 
those whose Nyquist rates were low. 

In Chnpter 3, Section 2 the coarse-sampled-data 
filter using a digital counter averager wat^ investigated. 

It was found that the output of this filter was equal to 
that or the linear digital filter in neither its expectation 

The coarse-samplal-^data filter 


nor its autocorrelation. 



wna n. citified to yield the N-bit feedback sarrilod- 
dfitri { lltf-r i jiveatiga feed in Oiwpter 3 , Section 3. It WiS 
fo'ml tlin t thtf expected value of thr^ output of the lottor 
* titer am that of the output of the linear digital filter 
weT"o equf'l. It was also found that the autocorrelation 
of the output of the N-bit feedback sampled-data filter 
approached asymptotically tho autocorrelation of the output 
of the linear digital filter as the dithering rate 
d/T' = K/T ) increase:, for a fixed inout sampling rate 
(l/T). One disadvantage of the N-Jbit feedback saniplcd- 
data filter compared to she coarse-sampled-data filter 
is that for a linear difference equation of order M, the 
corresponding coarse-sampled-data filter requires the 
length of the feedback register to be M, while its mod if i- 
oabldn requires the register length to be MN, where N is 
the mtio of the input sampling period T to the period 
of the dither T’. 

In Chapter 4, the results obtained by simulating 
first and second order N-bit feed bad arapled-data 
filters are discussed. It is seen that- within a certain 
amplitude range of the input, the Mit feedback samp! 
data filter responds as a linea-- filter. The filter 
requires a minimum 'run time' dv>ng which the nature! 
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roaponae due to the initial conditions is allowed to 
doony to a ne, ‘Eligible value, after which the power sp<-''tit. 
of till’ outputs of the N-bit feedback sampled-data filter nryj 
the linear digital filter may be compared'. It is seen th'*t 
for the nonlinear filter realised in hardware or simal'’led 
on the computer, there is a specific dithering rote for n 
fixed input sampling rate for which matching of spoctm 
occurs. Below this value of N, the number of points svem^rGi 
to obtain y(lT) is inadequate, and hence the autocorrelation 
of y(lT) is not equal to that of y(lT). Above this 
optimuio value of N, the error introduced due to truncation 
and/or rounding increases and causes the power spectrum 
y(lT) to increase above that of y(lT). 

Scope for Further Work : 

■ The ooarse-sampled-dn f '! filter proposed by Ki rlin 
has been modified to yield the N-bit sampled-data filter; 
•this uses a 1-bit quantizer. The results may be extended 
to the case when the quantizer has a greater number of 
quantization levels. A sensitivity analysis may be 
carried out to determine the perturbation in the power 
spectrum due to perturbation of coefficient values. 

A requirement of this nonlinear filter is tivat the 
input sampling rate be low. Hence the filter may bo used 
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In J^yatems that mnlti-rate with 
Invariant Multirate FiJ Bering requires that an analog? 
i liter ha ret^lnoed by a bank of multirate linear digital 
filters. The outp\it;B of the individual filters are aided 
after being appropriately delayed, to yield the o/p wnich 
is covariance invariant with the o/p of the analog filter* 
It would be interesting uo investigate the power spectrum 
of the CIMR filter when each of the digital filters 
in the bank, is replaced by an N~bit feedback sampled- 
data filter. 

As stated previously in Chapter 1 , simulators oj 
fading disperaiv-i channels require a great amount of 
digital hardware. The linear coi,.plex filters could be 

impTemented using N-bit feedback coaree-sampled-da ta 
filters. A simulation of such a nonlinear complex filter 
wou3 be highly informative of the nonlineai ilter s 


performance. 
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